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1. INTRODUCTION 

During SESAME 1979, thunderstorms were 
overflown (at 18 km) with NASA's WB-57F equipped 
with a three channel (visible, 6.7 micron, and 
11 micron) scanning radiometer, the Cloud Top 
Scanner (CTS) . In addition to the CTS data, cloud 
top height information is available from a down- 
ward pointing lidar on the same aircraft and from 
side views of the thunderstorms photographed every 
30 seconds from a secuud aircraft at 14 km. Also, 
for the case study day discussed here (June 7, 

1979), two geosynchronous satellites were operated 
in the short -interval (3 minute) mode in a synchron- 
ized fashion for the calculation of stereo heights. 

The general objective of the study is 
to examine and compare the high resolution aircraft 
observations (primarily the 11 micron channel) with 
GOES measurements to better understand and inter- 
pret the satellite observations of severe thunder- 
storms and to help determine sensing parameters for 
future satellite missions. Specifically, the 
objectives are to: 1) determine the relation be- 

tween height deviations and IR temperature devia- 
tions of overshooting thunderstorm tops; 2) 
determine the effect of spatial resolution on 
estimates of thunderstorm, height and on calcula- 
tions of thunderstorm top ascent rate; and 3) 
attempt to determine the relation of storm maximum 
height (as determined from satellite or aircraft) 
and cloud top ascent rates, to vertical velocities 
estimated from one-dimensional cloud model runs. 

2. MATURE STORM CHARACTERISTICS 

2 .1 General Characteristics 

Seven successive overpasses were 
flown over a mature thunderstorm with a ten minute 
return time between passes. Fig. 1 displays a 
summary of cloud top equivalent blackbody tempera- 
ture (Tgg) as viewed from the CTS (aircraft) and 
satellite with other pertinent data, including 
maximum radar reflectivity and height information 
from the aircraft lidar and satellite stereo 
calculations. Fig. 2 shows views of the storm at 
the time of the first overpass (line 3) from the 
CTS, satellite, and the Learjet aircraft at 14 km. 
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The storm began to build about 2140 
GMT (1540 CST) an hour before the first overflight 
(see Fig. 1). The aircraft overflights began 
about the time of maximum storm intensity as indi- 
cated by the radar and satellite IR data. The 
CTS data are represented in Fig. 1 by spatial 
averages over the coldest 1 km and 10 km squares. 

While the satellite is indicating a minimum T BB 
of 2 2 OK, significantly warmer, or below, the 
tropopause, the CTS IR data is clearly indicating 
Tgg's as low as 203K (-70C), with a 1 km average 
of 206K for line 3. These low values of T BB 
indicate thunderstorm tops significantly above 
the tropopause and this is verified by the air- 
craft lidar maximum height of 13.9 km and a 
satellite stereo estimate of 13.4 km. The average 
T BB over a 10 km square (approximately the satel- 
lite instantaneous field of view, IFOV) is much 
warmer and corresponds more closely to the satel- 
lite observation. The large satellite IFOV 
obviously leads to a significant underestimate of 
thunderstorm height . 
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2.2 


Description of Line 3 


2.3 


Cloud Model Analysis in Comparison to 
Aircraft and Satellite Observations 


Fig. 2 shows the various views of the 
storm at about the time of the line 3 overpass. 
The approximate lateral boundaries of the CTS and 
the position of the lower altitude aircraft (L) 
are displayed on the satellite images. A hand- 
remapped radar image is also shown. 


The CTS visible image shows over- 
shooting cumuliform features, each having its 
accompanying Tg B relative minimum, or cold feature. 
Individual turrets observed in the CTS image can 
be identified in the side view photograph. The 
CTS data and the side look photography are clearly 
showing a multicell thunderstorm, although the 
satellite IR image has only one Tgg minimum, because 
of its coarse resolution. The radar image also 
indicates only one rainfall core. The minimum 
t BB is 203K (-70C) at a lidar-determined height of 
13.9 km (see Fig. 3), which is UK colder than the 
tropopause and 5K colder than the ambient tempera- 
ture at that altitude. 


The downward pointing lidar gives 
cloud top height information along the centerline 
of the CTS image for comparison with the CTS IR 
data. Lidar data was only available for clouds 
above 12.8 km, i.e. above the tropopause. The 
lidar height-CTS T BB data are very well correla- 
ted (0.91) and have a regression slope of 9.2K 
km -1 , very close to adiabatic. The aircraft T BB - 
height relation is shown in Fig. 3 as the dash- 
dot line above the tropopause. At temperatures 
colder than -60C the CTS T BB absolute values are 
accurate to approximately +2K, so it is possible 
the cloud top temperatures are slightly colder than 
indicated here, Because the cumuliform turrets are 
very dense, the radiation emitted from the cloud 
top is coming primarily from the first 100 m. 

Both the slope of the relation and the deviation 
from ambient indicate little mixing with the 
stratospheric environment. This appears to agree 
with the aircraft overpass data presented by Pitts 
et al. (1975) and disagrees with the ideas of 
Roach (1967), who postulated intense mixing at the 
cloud edge and cloud top Tgg’s close to ambient. 



Flgu/ie 3. Local mounding with, one- dimensional D,HECT,0N 
model results plotted and presented In table InseAt. 


The aircraft overpass of line 3 was 
at the time of approximate maximum storm height. 

This is obvious from subsequent aircraft passes, 
although the satellite data shows minimum Tgg 
10-20 min later. A second updraft surge reached 
a lower maximum height at line 7 (Fig. 1), and a 
much weaker, separate storm (discussed in 
Section 4) reached about 9 km. A one- dimensional 
cloud model (Simpson and Wiggert, 1969) was run 
for various updraft radii (0.5 - 2.2 km) on the 
Dodge City, Kansas 2315 GMT sounding adjusted for 
local surface conditions in the area of the cloud 
(Fig. 3). The model run with an input updraft 
radius of 2.0 km approximately reproduces the 
observed maximum height (13.9 km) and minimum 
temperature (-70C) . The solid line to the right of 
the temperature profile in the troposphere is the 
temperature trace produced by the ascending parcel 
in the model cloud for the 2 km radius. Because 
of the effect of entrainment, model clouds with 
smaller updraft radii grow to lower maximum heights 
as seen by the dashed line (to the left of the 
ambient temperature) which connects the maximum 
height s /minimum temperatures from the model runs. 

The cloud parcel ascends nearly adiabatically through 
the neutral point to its maximum height. However, 
when comparing different storms reaching different 
^min ’ s » one must use a smaller lapse rate ('WK 
km“l) to convert to height variations. A knowledge 
of the actual cloud top T m ^ n can be converted to 
height by using the dashed curve, not the ambient 
trace, and not an adiabat. However, the satellite 
overestimates the T m i n and must be corrected. A 
crude correction based on the CTS/GOES comparison 
in this case would be 


T . * T - 15. (1) 

cloud sat 

From this estimate of T c i OU( j an estimate of the 
updraft intensity can be made using the cloud 
model runs. The table in Fig. 3 gives the model 
results and the maximum height information from 
three aircraft overpasses. The line .3 height 
and temperature is roughly met by the model run 
with an updraft radius of 2 km. The associated 
maximum vertical velocity (W max ) is 35 m s“^ 
at 275 mb. Therefore, although the storm is 
observed to penetrate only 1.5 km above the 
tropopause, it has, from this analysis, a very 
strong updraft, which agrees with the large 
radar reflectivities observed. The reason the 
cloud penetrates only the short distance above 
the tropopause is the effect of the relatively 
high static stability in the upper troposphere 
above 9 km which lessens the accelerating effect 
of the buoyancy term relative to water loading. 

The combined effect decelerates the model parcel 
significantly before hitting the tropopause. 
Therefore, in order to convert measures of storm 
maximum height to estimates of updraft intensity 
(a parameter better correlated with the occurrence 
of severe weather), the sounding structure must 
be taken into account beyond just a knowledge of 
tropopause temperature or height . 


2 .4 Satellite Observed Cloud Evolution 

The solid line in Fig, 1 (GOES minimum 
T bb ) is the satellite observed minimum Tgg over 
the area observed by the CTS. As the cloud 
passes its most intense stage, the coldest portion 
of the cloud shifts downwind over the central por- 
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tion of the anvil in the satellite image. This 
feature is indicated in Fig. 1 as the dotted line 
above the solid line from 2250-2310 GMT. From the 
side-looking photograph, the satellite cold point 
is related to smooth cirrus debris ^drifting' downwind 
at an altitude at, or just above, the height of the 
now subsided convective area. The cirrus is at, or 
just above, the tropopause. The cirrus is not a 
candidate for producing a downwind warm spot some- 
times observed in satellite images of intense 
thunderstorms because of its lack of height above 
the tropopause and the decrease with height of the 
ambient temperature above the tropopause (Fig. 3). 

From the side photographs and the CTS 
overpasses, it is obvious that the convective area 
continually regenerates on the southwest end of 
the storm. However, following the coldest portion 
of the cloud from the satellite observations, one 
might incorrectly infer updraft movement into the 
middle of the anvil. A second updraft surge is 
obvious in the CTS data (lines 7 and 8) and side 
photographs and can be detected as a separate cold 
point in the satellite IR data (2250-2330 GMT). 

3. GROWING STORM CHARACTERISTICS 

A sequence of passes was also flown 
over a growing, but relatively weak, cumulonimbus 
cluster that reached a minimum CTS Tgg of 235K 
(-38C) , or a height of about 9 km. The trend of 
minimum Tgg is shown in Fig. 4. Because of the 
rapidly changing features, continuity of individual 
towers is difficult in many cases (side photography 
was not available at this time) and Fig. 4 gives 
the minimum Tgg anywhere in the cluster at each time. 
The concurrent satellite observations indicate ob- 
served Tgg ’ s much warmer (30-40K) than that observed 
with the CTS. This is due to field of view (beam 
filling) effects, along with the satellite "seeing" 
the side of the cloud. 

The period of rapid cloud growth has 
a CTS-measured ascent rate of 2.2K min“l. The 
minimum Tgg (-38C) corresponds approximately with 
the cloud model run for 0.7 km radius (see Fig. 3). 
This model run has a vertical velocity peak of 
15 m s~l at 500 mb ('wlOC) . Over the temperature 
layer spanned by lines 10 and 11 (-16C to -38C) 
the average model w is 10.1 m s~l. Using the 
lapse rate following the parcel in the model run 



figure. 4 . Equivalent blackbody tempwoutuKe (Tgg) 
cu a function o& time glowing btoim a 6 obaeAved 
by aJjicACL^t and tateilite ladiomet&u * . 


(8.3K km"*)» the CTS ascent rate of 2.2K min“l 
converts to a w of 4.4 m s“^. Thus, the cloud 
top ascent rate appears to be only about half of 
the vertical velocity. This observation is in 
agreement with other studies for decelerating 
convective boundaries. 

4 . SUMMARY 

Comparison of thunderstorm top obser- 
vations made by aircraft overflights with a 
scanning radiometer and concurrent satellite 
observations indicate that: 

1. Overshooting, cumuliform features 
on various scales appear as distinct cold areas in 
the high resolution (80 m) IR (11 micron) aircraft 
images . 

1 . Primarily because of field of view 
effects, the satellite IR observations overestimate 
Tgg (minimum) by about 15K for storms significantly 
penetrating the tropopause. A satellite field of 
view of 1 km would alleviate this problem. 

3. Profiles of height and Tgg across 
overshooting features indicate an approximate adia- 
batic cloud surface lapse rate. However, one- 
dimensional cloud model runs point to using a 
lapse rate of ^7K km - ^ when comparing different 
Storms reaching various T m ^ n 's. 

4. In order to convert measures of 
storm maximum height to estimates of updraft 
intensity (a parameter better correlated with the 
occurrence of severe weather) the sounding structure 
must be taken into account beyond just a knowledge 
of tropopause temperature or height. 

5. During periods of updraft weakening, 
the lowest, satellite-observed Tgg may occur down- 
wind with cirrus debris. 

6. For small, growing storms, the 
satellite will overestimate cloud top Tgg by as 
much as 30-40K, again primarily because of field 
of view effects. 

7. The cloud top ascent rate of a 
growing storm is approximately half the magnitude 
of the vertical velocity calculated from a one- 
dimensional cloud model. 

5. FUTURE WORK 

A detailed description of aircraft-^ 
observed features and storm evolution will be 
completed for this case study. Also, techniques 
are being developed to correct satellite-observed 
^min* s f° r satellite field of view effects and 
other effects so that a better estimate of cloud 
top temperature (and height) can be made from the 
satellite data. The aircraft data, along with 
satellite stereo and high resolution IR data from 
low orbiting satellites will be used to verify and 
tune the techniques. 
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ABSTRACT 

Digital infrared data from a geosynchronous satellite (SMS 2) on 6 May 1975 are used to study thunder- 
storm vertical growth rates and cloud top structure in relation to the occurrence of severe weather (tor- 
nadoes, hail and high wind) on the ground. All thunderstorms from South Dakota to Texas along a north- 
south oriented cold front are monitored for a 4 h period with 5 min interval data. 

An examination of five cloud elements having eight tornadoes indicates that in seven of eight cases the 
first report of the tornado took place during, or just after, a period of cloud top ascent. This vertical velocity 
is applicable to an area of 1 5 km on a side. 

Thunderstorm growth rate, as determined by the rate of blackbody temperature isotherm expansion 
and minimum cloud top temperature, are shown to be correlated with reports of severe weather on the 
ground. A time analysis indicates that the derived parameters reach critical values soon enough to provide 
a potential warning lead time of approximately 30 min. 

Equations are derived relating the thunderstorm growth rate to vertical velocity and outflow layer 
divergence. Severe thunderstorm elements are shown to have mean vertical velocities approximately twice as 
large as the non-severe elements. The outflow layer divergence is calculated to be 1 X 10~ 3 s' 1 for the 
severe thunderstorms. 


1. Introduction 

Information obtained from geosynchronous satellite 
data has the potential of leading to greater under- 
standing of convective and other mesoscale activity 
and also the potential of being an important com- 
ponent in the detection and monitoring of thunderstorm 
activity in general, and of severe thunderstorms in 
particular. The detection of thunderstorms using satel- 
lite data has centered on the interpretation of images. 
For example, Purdom (1976) has noted that inter- 
secting cloud lines can produce enhanced convection 
and thunderstorms. However, quantitative information 
on cloud growth rates can also be determined by 
examining a sequence of images. Sikdar et al. (1970), 
Purdom (1971) and Am (1975) have used ATS 3 visible 
data to measure thunderstorm anvil expansion rates. 
Adler and Fenn (1976), Yuen (1977) and Negri et al. 
(1976) have used SMS data to make similar calcula- 
tions. Adler and Fenn (1976) showed examples of cold 
area expansion rates related to thunderstorm growth 
using SMS window channel infrared data. These 
measurements were for areas smaller than entire cirrus 


1 A preliminary version of this paper was presented at the 

Tenth Conference on Severe Local Storms at Omaha, Nebraska. 


anvils and provide a technique to monitor individual 
thunderstorms rather than a group of storms under 
a large cirrus shield. This last study also included 
an example of a tornado preceded by the rapid ex- 
pansion of a cold area, implying rapid ascent. All of 
these studies indicated that there are parameters de- 
terminable from geosynchronous satellite data that 
tend to be related to the occurrence of severe weather. 

2. Objective and approach of study 

The objective of the research reported on in this 
paper is to develop techniques to use geosynchronous 
satellite data to help detect and monitor severe 
thunderstorms, and to study their characteristics. 
Specifically, the purpose is to determine if thunderstorm 
growth rates and characteristics observable with SMS 
infrared (IR) data can be correlated with the occur- 
rence of severe weather on the ground. The occurrence 
of thunderstorm-related severe weather (tornadoes, 
hail, high winds) is highly correlated with the inten- 
sity of convection. Using SMS IR data one should be 
able to calculate parameters related to convection 
intensity, such as cloud top temperature and its rate 
of change. 

A particular area is outlined and all thunderstorms, 
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or thunderstorm clusters, in that area for a given 
time period are monitored. Parameters derived from 
the satellite data are then compared to ground reports 
of severe weather. The discussion of the results in 
this paper concentrates on examples of tornado- 
bearing thunderstorms and comparisons between severe 
and non-severe thunderstorms. Vertical velocity and 
divergence estimates are also presented. 

3. Data and method of calculation 

SMS 2 digital IR data from the period 1803-2208 
GMT 6 May 1975 were used in this study. This 
period was part of a limited-scan observation period 
wherein a relatively narrow swath, with its long axis 
oriented east-west, was observed by the satellite, but 
on a frequent basis. For this case the nominal period 
between observations was 5 min. The spatial resolu- 
tion of the IR channel is about 9.4 km at 40°N. 
There is a sampling overlap in the left-to-right direc- 
tion so that each data point is centered 4.7 km from 
the center of the next. Lines are repeated during the 
data processing to complement the left-to-right over- 
sampling so that image distances, either left-to-right 
or up-to-down, correspond approximately to each 
other in terms of geographical distance. Because of 
the sampling overlap in the left-to-right direction and 
the repetition of lines, each data point in this study 
represents an area approximately 4.7 km on a side, 
or 22 km 2 . 

The calculations presented in the next section were 
performed on the Atmospheric and Oceanic Informa- 
tion Processing System (AOIPS), an interactive image 
processing system described by Billingsley (1976). 
With this system, an irregularly shaped area can be 
outlined on a color enhanced television image and 
a histogram of the digital counts inside the area 
obtained. 

The temperature resolution of the IR data is 1 K 
in the range of equivalent blackbody temperatures 
(TWs) of interest in this study. Due to the fixed 
point algorithm used in the transformation of raw 
counts into calibrated values, voids appear in histo- 
grams (Goddard and Remondi, 1975). In the range 
of Tbb between 206-226 K, these voids appear every 
3-5 K. 

From the calculated histogram a cumulative histo- 
gram is compiled starting from the cold end of the 
temperature distribution. The number of points Ni 
with Tbb^Ti is obtained, with 7\ ranging from the 
lowest temperature in the area to the highest tem- 
perature To, where To is defined as the warmest Tub 
that forms a closed isotherm around the element 
being defined. In the case of easily definable, indi- 
vidual thunderstorms or thunderstorm clusters, an 
entire cirrus anvil may be defined. Examination of 
cirrus anvils led to the subjective determination that 
the 7^=258 K isotherm corresponds approximately 


to the anvil edge. Although this value is much warmer 
than the actual cloud temperature, the emissivity at 
the thin anvil edge is much less than unity, resulting 
in the warm T B b . In addition, the response function 
of the sensor may result in a more gradual change 
in Tbh than is actually occurring. Preliminary exami- 
nation of expanding anvils indicated that the areal 
change within somewhat lower Tbb 1 s paralleled the 
areal change within 258 K. For this reason a lower 
temperature, Tbb— 226 K, was chosen as the warmest 
cloud tops to be considered in the definition of areas. 
Therefore, the clouds studied here have already ob- 
tained a height of at least 9.5 km based on soundings 
for this day. The early development of the thunder- 
storms is thus not examined in this study. Satellite 
observations of earlier growth of thunderstorms is 
examined by Adler and Fenn (1978). 

As the anvil of a thunderstorm or thunderstorm 
cluster expands, it eventually comes into contact with 
other anvils. When the warmer temperatures near the 
outer edge of the cloud can no longer be uniquely 
identified with the feature being examined, T 0 is de- 
creased to a value such that the statistics of the data 
inside the area defined by T 0 now represent the feature 
being examined. In areas where a number of anvils 
have formed a continuous cirrus overcast, the object 
is to detect and monitor colder (higher) regions within 
the general cirrus background. In these cases To is 
often very cold (' — 212 K) and the elements defined 
contain only a small range of Tbb . 

The thunderstorms analyzed in this study lie along 
a generally north south oriented cold front from 
Nebraska southward through Texas. This area con- 
tains all the reported tornadoes for this day during 
this time, other severe thunderstorms, and thunder- 
storms with no reported severe weather. A selection 
of SMS 2 images covering the time period studied is 
shown in Fig. 1. Especially late in the period there 
are areas of cirrostratus clouds which cannot be asso- 
ciated with individual thunderstorm anvils. Thunder- 
storm elements can still be identified bv analyzing 
only the colder portions of a cloud. For example, 
Fig. 2 shows the distribution of selected Tbb isotherms 
for the region outlined in the 2128 GMT panel of 
Fig. 1. The area of interest in Fig. 2 is the region of 
relatively cold temperatures below the center of each 
panel. At 2057 GMT there is only a small area having 
a temperature of 208 K. Over the next 30 min the 
208 K isotherm expands rapidly and colder tempera- 
tures appear. These changes imply ascent. The element 
defined by the 208 and 206 K isotherms is typical of the 
small, cold cloud elements monitored in this case. This 
particular element (defined as cloud 8b) was associated 
with the Omaha tornado which touched down at 
2133 GMT. The position of the Omaha tornado rela- 
tive to the cloud top Tbb s is shown in the last panel 
of the figure. A hail report is shown in the panel for 
2118 GMT. A discussion of tornado reports relative 
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Fig. 1. Sequence of SMS 2 images over study area and period. 
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Fig. 2. Analyses of infrared T bb isotherms over a small area shown in Fig. 1. 
H marks location of hail reports; triangle marks tornado report. 


to cloud top features is given in the next section. 4. Results from tornado-bearing storms 
From Fig. 2 it can be seen that the expansion of the 

small cold areas cannot be interpreted as anvil ex- Five clouds or elements in the area examined had 
pansion. Analysis of the cold area expansions in terms tornadoes associated with them, with the total number 
of vertical velocity and divergence is given in Section 7. of tornadoes being nine. The approximate locations 
The occurrence of severe weather in this study is of the tornadoes are given in Fig. 3. Four of these 
determined primarily from the Severe Weather Events tornadoes can be considered relatively strong tor- 
Log compiled bv the National Severe Storm Forecast nadoes. These are identified by nearby town names 
Center (NSSFC) in Kansas City, with cross reference as the Pierce, Magnet, Winside and Omaha tornadoes, 
to the Environmental Data Service’s Storm Data. all occurring in Nebraska. The fifth tornado cloud is 
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Fig. 3. Location of tornadoes on 6 May 1975 in area of 
interest. Arrows represent approximate tracks for relatively 
strong tornadoes. 

associated with the relatively weak Saunders County 
(Nebraska) tornado. 

Fig. 4 is the thunderstorm growth rate diagram for 
cloud 4. The cloud numbering system is the original 
one used in the study to identify approximately 
60 clouds or elements from South Dakota to Texas. 
Eventually, only 39 elements were determined to be 
active thunderstorm areas. The eliminated elements 
could only be followed for a short time and indicated 
little or no growth. The cloud or element numbers 
will, therefore, not necessarily be sequential. The 
diagram indicates the relation between N t and time, 
where is the number of data points in the defined 
element with blackbody temperature Tsb^T^ For 
example, the curve in Fig. 4 for 7\ = 218 K exhibits 
an increase from .V=14 at 1808 GMT to iV=60 at 
1823 GMT. Before 1808 the minimum Tub was 
warmer than 218 K. The curve is drawn to N= 1 at 
1803 GMT, so that the growth from zero to 14 can 
be displayed. Because the A r axis is logarithmic, the 
slopes of the lines are proportional to N~ l dN/dt. 
Plotted points which are circled indicate that the 
area is not defined uniquely by that isotherm, but 
that a portion of the border is approximated. The 
approximation is usually small and does not signifi- 
cantly affect the count. The conversion from number 
to area is accomplished by multiplying the number 
of points by 22 km 2 . The location in time of the 


severe weather events is shown across the top of 
the diagram. 

Cloud 4 in Fig. 4 shows a rapid decrease in mini- 
mum blackbody temperature, or T min , and a rapid 
expansion of areas within isotherms. The decrease in 
temperature can be calculated by moving horizontally 
across the diagram at, for example, iV=10. The 216 
and 212 K lines are separated by 20 min, indicating 
a rate of decrease of 0.2 K min -1 . The values of 
N^dN/dt are of order of magnitude 10 -3 s' 1 . The de- 
crease in temperature and expansion of isotherms, of 
course, implies ascent. The vertical velocity is related 
to the rate of temperature change by 



With a reasonable lapse rate, 8 K km -1 , the tempera- 
ture change of 0.2 K min' 1 is equivalent to a vertical 
velocity of w of 0.4 m s _1 . This is a reasonable vertical 
velocity estimate considering the area over which it 
is an average — 225 km 2 ) and considering it represents 
a time when the cloud top is decelerating as it ap- 
proaches and penetrates the tropopause at approxi- 
mately 12 km. Much larger values of vertical velocity 
are calculated for clouds whose tops are just 1 or 2 km 
lower. A more detailed examination of the vertical 
velocity estimates is given in Section 7. 

Three tornadoes are associated with cloud 4 (Fig. 4). 
One in Bon Homme County (South Dakota) is in the 



TIME (GMT) 

Fig. 4. Thunderstorm growth rate diagram for cloud 4. N is 
the number of data points in the element with blackbody tem- 
perature TsB^Ti, where 7\ is as identified in the diagram. 
Occurrence of severe weather is marked by arrow at time of 
report, with T indicating tornado. 
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Fig. 5. Thunderstorm growth rate diagram for cloud 5. 


vicinity of this element, but when no distinct cold 
area can be defined. The Pierce tornado, the strongest 
of the three, was first reported at 1905 GMT. At that 
time and immediately preceding it, during the tornado 
formation period, the element exhibits cold area, ex- 
pansion and decreasing temperature, indicating ascent. 
The same pattern also is associated with the Knox 
County tornado. 

Fig. 5 shows the diagram for cloud 5 which is 
identified as having been associated with the Magnet 


tornado. At the beginning of the period, relatively 
slow growth is indicated for the 221 and 220 K areas. 
A few points with T=218K appear and then dis- 
appear before a period of very rapid growth is ob- 
served between 1825-1840 GMT. The maximum of 
N-'dN/dt is 6.2 X10 -3 s _1 at 218 K. The minimum 
temperature also drops very rapidly during this period. 
After 1843 GMT no clearly defined cold area can be 
identified, although the cloud element can still be 
identified in the higher resolution visible images. In 
the IR data cloud 5 appears as an extension of an 
element immediately to the south, cloud 22. Not until 
1923 GMT or 40 min later, is a definable cold area 
at 211 K observed for cloud 5. Over the next 10 min 
this isotherm expands swiftly. The initial report of 
the tornado is at 1915 GMT, as noted on the diagram, 
but the tornado hits the town of Magnet at 1945 GMT 
and lasts at least until 2000 GMT. It, therefore, ap- 
pears that this tornado also formed and touched down 
during a time when the thunderstorm top was as- 
cending at the spatial scale observable with the 
IR data. 

Just to the south of the Magnet cloud is cloud 22 
(Fig. 6). This element can be defined continuously for 
almost 3 h. Two tornadoes are associated with this 
element at two widely separated times. The first, at 
Winside, Nebraska, is reported at 1945 GMT. As with 
the tornadoes already examined, this tornado occurs 
during a period of rapid increase in cold area, which 
in this case lasts from about 1900 GMT until 2000 
GMT. The maximum value of N~ l dN/dl is 4.5 X 10~ 3 s -1 , 
following the 211 K isotherm. The second, weaker 
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Fig. 7. Thunderstorm growth rate diagram for cloud 38. 

H indicates report of hail. 

tornado, at 2i30 GMT, is located in Dixon County, 
Nebraska. Although there is a brief period of expanding 
cold area between 2125-2145 GMT, this tornado 
occurs during a period of generally decreasing cold 
areas and warming temperatures. One can speculate 
that during the dying phase of element 22, one last 
vigorous thunderstorm penetrates the cirrus shield, 
showing up for a short time as the secondary maxi- 
mum in the 209 K curve at 2140 GMT. 

A tornado occurring in Saunders Country, Nebraska, 
is associated with cloud 38 (Fig. 7). Again, the reported 
time of the tornado is during a period of cold area 
expansion indicating ascent. The most rapid increase 
in N is between 1930-1950 GMT at a rate of 6.7 
X 10 -3 s” 1 . The minimum temperature reached is 
208 K. 

The curves for cloud 8b are shown in Fig. 8. Two 
tornadoes are associated with this cloud, including 
the severe Omaha tornado. This element is observed 
to go through two periods of rapid ascent during 
a 2 h period. Reports of hail and high wind (>50 kt) 
are associated with the first period of r&pid expansion. 
The element is then observed to have a decrease in 
cold area between 2040-2055 GMT. After that there 
is a sharp increase in A (208 K) and colder tempera- 
tures appear. Although there is good continuity of 
the element during this cycle of inferred ascent and 
descent, the second period of ascent is probably 
brought about by a new thunderstorm penetrating 
the cirrus shield in approximately the same area. 

Similar to previous cases, the growth rate of cloud 8b 


during the 30 min period prior to the first tornado 
report indicates ascent. During the lifetime of the 
tornado (2133-2150 GMT) no very sharp area in- 
crease is noted. In fact, the dip in the 206 K line 
could possibly be interpreted as a partial collapse of 
the thunderstorm top. This is the only element for 
which such a feature is evident. No data are available 
from 2150 to 2208 GMT. At the later time, 205 K 
temperatures appear, indicating more growth. The 
Beebeetown, Iowa, tornado is reported during the 
period 2200 to 2230 GMT. 

During the time period examined, nine tornadoes 
were reported, eight of which could be clearly con- 
nected with a small, colcl area in the IR data. The 
exception was the Bon Hbttrme County tornado which 
was located in the vicinity of cloud 4, but long after 
it was no longer possible to define a small, distinct 
cold area in it. Of the remaining eight cases, all but one 
(Dixon County — cloud 22) occurred during, or just 
after, a rapid increase in cold area, indicating cloud 
top ascent. Therefore, it appears that the formation 
of tornadoes in a large majority of cases occurred 
during a period of increasing thunderstorm cloud top 
height. This upward vertical velocity applies on a 
spatial scale of approximately 15 km on a side. 

The observations of increasing thunderstorm top 
heights at the time of, or just prior to, tornado 
touchdown in the present study seemingly contradict 
aircraft observations of overshooting tops discussed 
by Fujita (1973), Pearl (1974) and Umenhofer (1975). 
They relate collapsing cloud domes (above the anvil 
cirrus) to the occurrence of tornadoes. This discrep- 



Fig. 8. Thunderstorm growth rate diagram for cloud 8b. 
W indicates high winds. 
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Fig. 9. Location of tornado reports relative to IR cloud top features. 


ancy may be related to a difference in scale between 
the two types of observations. The satellite observa- 
tions apply to a scale of about 15 km on a side. Shenk 
(1974) has examined the size distribution of over- 
shooting domes or turrets and finds a median diameter 
of approximately 5 km. Therefore, the satellite vertical 
motions are applicable to an area approximate!)' an 
order of magnitude larger than that observed with 
the aircraft. Comparison of results is therefore difficult. 

Radar observations also do not show a sharp in- 
crease in echo top height prior to tornado touchdown. 
For example, Lemon et al. (1978) show nearly steady 
or slightly decreasing radar heights for the 30 min 
prior to the Union City, Oklahoma, tornado. A rapid 
decrease of the radar top follows the tornado touchdown. 
Again these observations are on a small scale com- 
pared to the satellite-based cloud top changes. A full 
understanding of the relation of satellite-inferred cloud 
top changes, radar echo height and actual cloud 
changes on various scales must await closer coordina- 
tion in the observations. 

Locations of six tornado reports relative to cloud 
top features observed in the satellite IR data are 
shown in Fig. 9. A similar diagram for the Omaha 
tornado is given in Fig. 2. The surface tornado loca- 
tions have been shifted to cloud top relative positions 
for a cloud height of 12 km to compensate for the 


satellite viewing angle. The tornado reports tend to 
be located to the west or southwest of the coldest 
temperatures, in an area of large gradient of T bb> 
Early in the 4 h period (before 2000 GMT) the tornado 
reports are located near the 226 K Tim isotherm, 
which is near, but inside, the upwind cirrus anvil 
edge. Later, after the cloud has expanded and the 
anvil edge is farther from the cold center, the tornado 
reports are located more distant from the cloud edge. 

5. Example comparison of a severe and non-severe 
case 

Two adjacent thunderstorm elements which de- 
veloped in southeast Nebraska at about the same 
time provide an excellent comparison of clouds with 
and without accompanying severe weather. The dia- 
grams for these elements (clouds 8 and 9) are shown 
in Figs. 10 and 11, respectively. 

Both elements reach a temperature of 226 K a little 
after 1800 GMT and continue to grow to higher 
heights and colder temperatures during the next two 
hours. However, there are significant differences in 
the rate of growth and the final temperature reached. 
Cloud 9, with no accompanying severe weather reports, 
has lower values of normalized growth rate N~ l dN/dl 
and an associated less rapid decrease of minimum 
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temperature. Its eventual coldest temperature is 212 K 
compared to 209 K for cloud 8. Two reports of hail 
and two reports of high wind are associated with 
cloud 8. The first pair of hail and wind reports occur 
near the end of the period of rapid growth ; the second 
pair happen during a relatively stable period. This 
second pair may be due to a thunderstorm growing 
rapidly from below, but hidden by the already high 
and dense cirrus anvil of the previous cell. Finally, 
at the end of the period the new thunderstorm pene- 
trates the anvil and becomes observable as the 209 K 
area at 2015 GMT. 

The comparison of cloud elements 8 and 9 is the 
best available of those analyzed. These two clouds 
began about the same time, lasted a substantial 
amount of time and were easily defined. Other ele- 
ments were sometimes obscured by anvils of more 
rapidly growing storms and other weak elements often 
became undefined in a general cirrostratus shield. 
In general, however, the strong thunderstorms are 
easily observable. 

6. Results for all storms 

Results from the analysis of 39 thunderstorms or 
thunderstorm clusters are discussed in this section. 
There are 15 elements associated with severe weather 
and 24 with no reported severe weather. In the area 
analyzed these 39 elements are all the clouds defined 
as growing thunderstorms. A few other elements were 


defined but lasted only a short time, showed minimal 
change, and no severe weather reports were associated 
with them. 

An examination of the thunderstorm growth dia- 
grams and the patterns of digital data gave the im- 
pression that two simple characteristics were correlated 
with the occurrence of severe weather. These charac- 
teristics are minimum cloud top temperature and 
thunderstorm growth rate as given by the normalized 
rate of cold area expansion. 

Fig. 12 displays the results of plotting the minimum 
temperature during an element’s lifetime (T min ) against 
the thunderstorm growth rate. The dots represent 
non-severe cases, while T’s and H’s represent tornadoes 
and hail respectively. The following procedure was 
used to arrive at a single value of N~ l dN/dt to repre- 
sent each element. The maximum N~ l dN/dt value, 
for any temperature ^226 K, is determined by ex- 
amining the time change between A 7 — 6 and /V = 30 
for each temperature. These two values were chosen 
so that the rate of increase is for relatively small 
areas, but that there are enough data points (6) to 
establish the lower threshold without being affected 
severely by instrument noise. A few elements had no 
curves of N lasting from N = 6 to 30, so the parameter 
was calculated using the largest range available. 

An examination of Fig. 12 indicates that the severe 
weather elements tend to have cold minimum cloud 
top temperatures and large rates of growth. A series 
of discriminant analyses (Panofsky and Brier, 1963) 
were performed, using the data displayed in Fig. 12. 
Table 1 shows the results. For 7\ nin by itself a severe, 



Fig. 11. Thunderstorm growth rate diagram for cloud 9. 
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non-severe dividing line of 212.8 K was derived. For 
temperatures above 213 K only one out of 12 elements 
had severe weather, while 14 out of 27 with tem- 
peratures colder than 213 K had associated severe 
weather reports. For maximum rate of growth the 
derived dividing line is 3.5X10 -3 s _1 with severe and 
non-severe events distributed around that value as 
shown in the table. A combined discriminant analysis 
using T m m and N~ l dN/dt produced the results shown 
in the last row in the table and by the diagonal line 
in Fig. 12. Based on the two variable discriminant 
analysis a probability of detection (POD) equal to 0.73 
was calculated along with a value of 0.31 for the 
false alarm rate (FAR) and 0.55 for the critical success 
index (CSI). The definition of these variables follows 
Donaldson et al. (1975). There are five elements with 
associated tornadoes, and their relations to the derived 
critical values are given by the fractions in parentheses 
in Table 1. 

Statistical tests were also performed to determine 
if there was a significant difference between the means 
for severe thunderstorms and non-severe thunder- 
storms. For T min the mean values were 209.5 and 
212.3 K for severe and non-severe cases. Using a 
/-distribution test, the difference of means was shown 
to be significant at the 5% level, with the value of / 
just missing the value for the 1% level. For N~ l dN/dt 
the means were 4.5X10' 3 and 3.1X10 -3 s~ l . This dif- 
ference of means also is significant at the 5% level, 
again barely missing the 1% value. 

It is not surprising that F min and thunderstorm 


Table 1. Discriminant analysis summary for cloud top mini- 
mum temperature, maximum rate of growth, and the two variables 
combined. The parameter L is positive for predicted severe 
weather. 




Fraction of elements 
with severe weather* 

Variable (s) 

Value for L =0 

L <0 

L >0 

Cloud top minimum 
temperature (Tmin) 

212.8 K 

12 \12/ 

27 \n) 

Maximum rate of 
growth 

(±4Z.\ 
\N dt ) 

3.5X10-* s-» 

22 \22/ 

-f-') 

17 Vl7/ 

Combined T m in and 
1 dN 
N dt 

*-•»+*« Gif) 

—0.048 Tmln 

in 

23 \23 / 

iifiA 

16 Vl«/ 


* Values in parentheses indicate tornado cases. 


growth rate ( N~ l dN/dt ) are correlated with occur- 
rences of severe weather. Both parameters are obvious 
indicators of the intensity of convection. Colder tem- 
peratures imply higher thunderstorm heights, which 
have long been related to thunderstorm severity 
through radar observations. All 39 elements involved 
in this analysis were associated with vigorous thunder- 
storms. The very weakest, as indicated in the bottom, 
left-hand corner of Fig. 12, had a minimum tempera- 
ture of 219 K, which is equivalent to a height of 
11.1 km, using the nearby 0000 GMT 7 May Omaha 
sounding. The coldest cloud had a temperature of 
205 K, or a height of 13.0 km. The 0000 GMT 7 May 



Fig. 12. Relation of thunderstorm rate of growth (cold area expansion) to cloud 
top minimum temperature and the occurrence of severe weather for all 39 thunder- 
storm elements in analysis area. T indicates tornado, H represents reported hail. 
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Table 2. Median warning lead time for all 
elements with severe weather. 




Median value of 


Number of 

lead time 

Criterion 

cases 

(min) 

Tmin 

14 

24 

1 dN 
N dt 

10 

27.5 

Combined 

11 

7 


sounding has a minimum temperature of 206 K. The 
tropopause is probably located at 222 mb at a tem- 
perature of 215 K, although there is another sharp 
increase in static stability above that, at a pressure 
of 168 mb and at a minimum temperature of 206 K. 
Twelve hours earlier the Omaha sounding has the 
tropopause located at 200 mb and 213 K. Farther 
south along the area studied the tropopause tempera- 
ture was also in the vicinity of 213 K. At Monett, 
Missouri, the tropopause temperature values were 211 
and 215 K at 1800 GMT 6 May and 0000 GMT 
7 May. Thus an approximate tropopause temperature 
for this region on this day was 213 K, which is the 
threshold temperature in the discriminant analysis 
discussed earlier. This result may be accidental, but 
is supported by similar observations by Fry or (1978). 
Additional studies are needed to ascertain the rela- 
tions among the temperature structure and tropopause 
height, the satellite-observed cloud-top temperatures 
and storm severity. 

The thunderstorm heights estimated from the satel- 
lite data in this study are substantially lower than 
the corresponding radar estimates. A cursory com- 
parison of the satellite-based heights with the radar- 
estimated storm tops as analyzed on the National 
Weather Service’s radar summary chart, indicates 
a mean 2 km underestimate from the satellite 
information. 

The reason for the radar-satellite differences is 
probably related to the scale of the observations and 
the limitations of the satellite data. The satellite- 
observed cloud top temperatures are for a larger 
horizontal area. The instantaneous field of view 
(IFOV) of the satellite IR sensor at this latitude is 
approximately 9 km on a side. The radar observations 
are probably representative of a smaller region. 

The satellite data are also limited by the tempera- 
ture structure. For elements of the size with which 
we are dealing, the cloud top probably takes on a 
temperature which is a combination of updraft tem- 
perature (from the moist adiabat) and the ambient 
air at that altitude. This will have the effect of 
somewhat limiting the greatest ideated height. 
Clouds penetrating into the lower Stratosphere will 
have their height underestimated. 

In addition, Negri el al. (1976) have already noted 


that SMS IR temperatures underestimate thunder- 
storm height, especially for small elements. They 
attribute this effect to inadequate sensor response 
when going from a warm to a cold target. This is 
also a probable contributor to the radar-satellite dif- 
ference noted in this study. Because of these factors 
the satellite-based estimates of thunderstorm height 
should be treated cautiously, although relative heights 
are probably valid. 

The thunderstorm growth rate parameter N~ l dN/dl 
is an indicator of the magnitude of the upward vertical 
velocity and the upper level divergence. The relation 
of these three variables is explored in the next section. 
One interesting point from Fig. 12 is that r min and 
N~ l dN/dl appear to be almost uncorrelated. One 
would expect that faster growing storms would pene- 
trate to higher heights and colder temperatures, and 
therefore that the two parameters would be correlated. 
It is believed that this lack of correlation arises 
because the thunderstorm elements in many cases 
were hidden by thick cirrus produced by other con- 
vection and, consequently, could not be monitored 
early in their history when the growth rates were 
most likely the most rapid. By the time the element 
penetrates the cirrus shield and growth rates can be 
calculated, the rate of growth has slowed. Such ele- 
ments reaching cold temperatures, however, are repre- 
sen Led by points on Fig. 12 with low temperatures 
and low growth rates. Another problem is the com- 
parison of N~ l dN/dl values irrespective of at what 
temperature or height they are calculated. Thus 
values of N~ l dN/dt of 5X10” 3 s" 1 calculated at 
Tub — 220 K may not reflect as large a relative velocity 
as 5X10 -3 s“ l at 210 K. With the addition of more 
cases these types of refined analysis will be possible. 

The parameters calculated appear to have the 
potential to be positive contributors to a severe storm 
warning system. To test what type of warning lead 
times might be achieved, the following calculations 
were made. Based on the criteria in Table 1, the time 
difference between when the element met the criterion 
and the time of the first report of severe weather was 
calculated for each element that met the criterion. 
The median values are displayed in Table 2. The 
median lead times are ~25 min for the one-variable 
criteria, but only 7 min for the combined criterion. 


Table 3. Median warning lead time for 
elements with tornadoes. 


Criterion 

Number of 
cases 

Median value of 
lead time 
(min) 

Tmin 

5 

33 

1 dN 



4 

35 

N dt 

Combined 

5 

29 
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A second calculation was made using the first report 
of tornadoes only and the results are shown in Table 3. 
The median lead time is ~30 min. The results of 
Tables 2 and 3 point out that the parameters cal- 
culable from the SMS data are not only correlated 
with reports of severe weather, but usually meet the 
severe thunderstorm criteria before the actual severe 
weather (hail, tornadoes) occurs. 

7. Vertical velocity and divergence estimates 

a. Interpretation of N~ l dN/dt 

The parameter chosen in this study to indicate 
thunderstorm rate of growth is the areal expansion 
of blackbody temperature isotherms in the window 
channel infrared data. This parameter has been chosen 
in place of rate of temperature change because many 
instances arise where this rate is difficult to calculate 
due to poor temperature resolution. Figs. 7 and 8 are 
examples of this situation. The calculated values of 
N~ x dN/dt also appeared to be relatively independent 
of the temperature or height. 

If a blackbody temperature isotherm is chosen so 
that it nearly coincides with the edge of the thunder- 
storm anvil, normalized expansion of the area within 
the isotherm is a direct measure of outflow divergence. 
That is, 

1 dA 

D (divergence) = , (2) 

A dt 

disregarding any dissipation of the anvil edge. In the 
present study, however, the isotherm expansions 
( N~ l dN/dt ) calculated were usually chosen for areas 
well within the anvil edge. There are two reasons for 
this choice. First, the emphasis here is on obtaining 
parameters which are applicable on a spatial scale as 
close as possible to thunderstorm scale. This is why the 
N~ l dN/dt values in Fig. 12 are calculated between 
N= 6 and iV =30. The anvil would cover an area 
two orders of magnitude larger. Second, it is often 
difficult to define an anvil edge because of complex 
interactions between thunderstorms. 

The calculations of N~ x dN/dt from N = 6 to A = 30 
are applicable to an area of about 15 km on a side 
and definitely do not represent anvil expansion. The 
results of this section will show that vertical velocity 
estimates can be made from the isotherm expansion 
values. 

The interpretation of N~ l dN/dt in terms of vertical 
velocity and divergence is shown schematically in 
Fig. 13. The cross-sectional or side view in upper 
part of the diagram shows the curved surface of the 
cloud top at two times, h and t 2 . The vertical co- 
ordinate is temperature. The emissivity of the cloud 
top is assumed to be unity. This should be a good 
assumption in the inner part of the thunderstorm top 
where this analysis applies. From times t\ to h the 



Fig. 13. Schematic views of idealized cloud top at two dif- 
ferent times. Temperatures are indicated by To, T\ and TV 
Solid arrows represent actual air motion, while dashed vectors 
represent the expansion of blackbody temperature isotherms, 
as viewed by the satellite. 

area inside the T 0 isotherm expands. This is seen by 
the horizontal dashed arrows in the upper part of 
the figure and is also indicated in the top view in the 
lower part of the diagram. This areal expansion is 
equivalent to the N~ l dN/dt which is measured with 
the satellite data. 

The actual flow of air, however, is given by the 
solid arrows. The air ascends and therefore cools 
(moist adiabatically) and is somewhat divergent. If 
this simple schematic diagram is basically correct then 
the expansion of cold areas as denoted by N~~ l dN/dt 
is really a combination of vertical velocity and di- 
vergence effects. 

b. Derivation of equations 

The relationships of N~ l dN/dt , divergence and 
vertical motion can be shown through the following 
derivation. A point on the cloud surface is represented 
by the position vector 

r=^i+yj+zk, (3) 

and the motion of the point is given by 

dr dx dy dz 

— = V = — H — jH — k. (4) 

dt dt dt dt 

The horizontal motion in the x direction (Fig. 13a is 
in the x-z plane) of a point on the cloud surface is 
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given by 


dx 
u — — 
dt 


dx 

dx 

— 

+w — 

dt 

* dz 


( 5 ) 


where dx/ dt | * is the horizontal movement of the cloud 
edge at a particular altitude or temperature. In Fig. 13a 
this is represented by the horizontal, dashed vector. 
The term dx/dz \ t is the slope of the cloud surface. 
Differentiating Eq. (5) with respect to x and adding 
the result to a similar expression for v differentiated 
with respect to y results in 


is derived. Eq. (10) states the relation among di- 
vergence, isotherm expansion, and vertical velocity. 
As will be shown, divergence is small compared to 
the other terms. The parameter s is large when the 
cloud has steep slopes and small when the slope is 
shallow. An expression for s, which can be used with 
SMS/GOES IR data, is derived in the next section. 
Although not used in this study, information from 
the higher resolution visible channel could be used 
to help infer cloud top shape during daylight hours, 
especially near sunrise and sunset. 


du 

dv 

u 

f dx > 

\+-l 

<dy \ 

dx 

dy 

dx > 

<dt z / 

' dy' 

\dt J 


dr 


+wV 2 — , 
dz 


( 6 ) 


if w is assumed to be constant in the horizontal. 
The left-hand side of Eq. (6) is the horizontal di- 
vergence D and the sum of the first and second terms 
on the right side is equivalent to the rate of isotherm 
expansion and the sum is equal to N~'dN/dt . Therefore 
Eq. (6) can be expressed as 


c. Formulation for parameter s 

Information concerning the shape of the cloud top 
and therefore the shape parameter s can be determined 
from the satellite IR data. The term in Eq. (9) is 
first expanded to obtain 

dr d /dx I \ d 

V 2 = — (— )+— 

dT dx\dT \ t / dy 



1 dN dr 

D = +w\ V— . (7) 

N dt dz 

The second term on the right side of Eq. (7) is 
a product of the vertical velocity and a divergence 
term which is a function of the slope of the cloud 
surface. By incorporating a vertical lapse rate into 
the expression, the vertical coordinate is shifted to 
temperature in Eq. (8) : 


Assuming symmetry in the two directions results in 
dr d /dx 


dT 


=H- \ 

dx\dT J 


The slope dx/dT is approximated by 
dx £N W * 
dT 2AT ’ 


(13) 


(14) 


1 dN / dT\ dr 

Z>= id )V 2 . (8) 

N dt V dz) dT 

Information about the shape of the cloud top can be 
inferred from the satellite data, by determining the 
area encompassed by the then coldest isotherm when 
the next, still colder isotherm appears. If the area 
of the first isotherm is large the slope is small and if 
the area is small, this implies a relatively steep slope. 
By defining a shape parameter j, such that 


where £ is the length represented by each satellite 
data point or pixel and N w the number of data points 
at T i when the next coldest temperature 7" n first 
appears. The expression assumes a square area. There- 
fore £NJ/2 is an estimate of the distance from the 
center of the isotherm area to the edge of the area. 
Divided by AT this produces an estimate of the slope 
in Eq. (14). The term AT is not always 1 K but is 
often approximated by 1.2 or 1.3 K because of the 
temperature void problem mentioned in Section 3. 

The finite difference approximation for Eq. (13) 
is simply 


r dr t 1 



v 2 — , 

(9) 

d /dx \ 

L arJ 

2 — (— ) 



dx\dT/ 


(dx/dT)!- (dx/dT) 2 


(15) 


one modifies Eq. (8) to become 


1 dN 

D= 

N dt 



( 10 ) 


Using Eq. (1) to substitute into Eq. (14) a useful 
alternate expression, 


1 dN 1 dT 

D— 1 , 

N dt s dt 


where the derivatives at positions 1 and 2 are the 
slopes at opposite sides of the cloud top, and Ax is 
the width of the area enclosed by the isotherm. As- 
suming they are equal in magnitude and opposite in 
sign and substituting Eq. (14) for the slope gives 

d /dx I \ 2 £NJ 

2 — ( — . (16) 

dx\dT \ t / AT Ax 

Because the calculations to follow are for areas of 


( 11 ) 
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the size of N= 6 to A = 30 satellite data points, an 
intermediate value of N — 16 was chosen as an ap- 
propriate size for the application of these formulas. 
Therefore, Ax in Eq. (16) is assigned a value of Vl6£ 
or 4{. Thus 

dr d /dx\ Nj 



A number of assumptions and approximations were 
needed to arrive at Eq. (17). A preliminary set of 
calculations using Eq. (11) with s equal to the inverse 
of NJ/2AT from Eq. (17) resulted in an average 
value of divergence D unrealistically small (5X 10~ 5 s _1 ). 
The calculations also indicated that because D was 
the small difference between the two relatively large 
values [see Eq. (11)], a small bias in s could drasti- 
cally change the results for divergence. Therefore, the 
parameter s is defined as 

2AT 

5 = a , (18) 

Nj 

where a is a constant to be determined empirically 
through an analysis of Eq. (11). 

The value for a was determined in the following 
manner. The terms N~ l dN/dt and dT/dt in Eq. (11) 
could both be calculated from the satellite data in 
a limited number of cases. Thirty-eight sets of vari- 
ables were available for 19 elements. The parameter s 
could also be calculated for these 38 cases from 
Eq. (18), at least in terms of a. The average diver- 
gence D in Eq. (11) was calculated using the con- 
tinuity equation for an incompressible fluid, 


(flatter tops) predominating at lower temperatures 
(higher heights), as expected. The fitted curve was 
determined by first dividing the data set into five 
approximately interval-equal groups based on tem- 
perature. Average s and temperature were calculated 
for each temperature interval and the linear regression 
curve was based on these five average points. This 
regression approach was adopted to reduce the effect 
of the relatively large number of points in the center 
region of the temperature range. 

d. Calculation of vertical velocity 

In Eq. (10), divergence is small compared to the 
other two terms, and to a rough approximation can 
be ignored to allow a direct relation between w and 
N~ l dN/dt. However, the effects of divergence in 
Eq. (10) can be partially taken into account by 
noting that, in the 38 cases previously examined where 
we estimate all of the terms in Eq. (10), D and 
N~ l dN/dt appeared to be correlated. This is reason- 
able since D is the horizontal expansion rate of a 
physical area and N~ l dN/dl is the expansion rate of 
area within Tim isotherms. Thus D can be approxi- 
mated by 

1 dN 

, ( 22 ) 

N dt 

where (3 is set at 0.15, based on the average values 

for D and N~ l dN/dt for the 38 cases. Substituting 

Eq. (22) into Eq. (10) and solving for w results in 

0.85s 1 dN 

w — . (23) 

( — dT/dz) N dt 


drt 

D = , ' (19) 

dz 

and a composite w profile derived from the 38 cases 
using Eq. (1). The resulting D is 0.5X10“ 3 s“ l . The 
mean value for N~ l dN/dt over the 38 cases is 3.4 
X10~ 3 s _1 . For s~ l dT / dt the average magnitude is 
(3.4Xl0“ 3 )/<*, where a is the proportionality constant 
from Eq. (18). Using these average values to substitute 
into Eq. (11) results in 

3.4 

0.5 = 3.4 , (20) 


3.4 

a = — = 1 . 2 . ( 21 ) 

2.9 

With a set equal to 1.2, the shape parameter s can 
be calculated from satellite data from Eq. (18) in 
some cases. Fig. 14 displays the distribution of 5 as a 
function of temperature for the 38 cases. The values 
span an order of magnitude with smaller values 


Eq. (23) provides a quick method of estimating 
vertical velocity even where there are no rapid de- 
creases in temperature, or where an element can only 
be viewed for a short time. 

Calculations of w using Eq. (23) are made for each 
of the 39 elements at particular temperatures. A lapse 
rate midway between a smoothed ambient sounding 
and a constant (9 K km -1 ) moist adiabatic rate is 
used. The 5 values are selected from the fitted curve 
in Fig. 14 for the various temperatures. For tempera- 
tures <T min a w of zero is assigned. The calculated 
w values are then averaged for the severe and non- 
severe elements and displayed in Fig. 15. The curves 
are fitted subjectively. 

The mean w for the severe elements is over 2.3 m s“* 
at 226 K or 10 km and decreases rapidly with height. 
The non-severe elements have nearly a factor of 2 
smaller vertical velocity values. Divergence values 
calculated for the two curves from 10-12 km using 
the continuity equation are 1.0X10 -3 and 5.0X10" 4 s -1 
for the severe and non-severe elements, respectively. 

These vertical velocity and divergence estimates are 
reasonable when it is recognized they are applicable 
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Fig. 14. Cloud top shape parameter s as calculated using 
Eq. (18) as a function of temperature. The 38 points are from 
19 elements. 

to areas ^15 km in a side. Vertical velocity estimates 
are very scale-dependent, and it is difficult to find 
values applicable on a scale equal to that determined 
from the satellite data. Cloud top vertical velocities 
estimated from aircraft observations and radar data 
provide at least a starting point. Observing penetrating 
or overshooting domes from aircraft produces vertical 
velocity estimates of the order of 1 m s -1 . For example, 
Umenhofer (1975) calculates the asceftt of an over- 
shooting dome to be approximately 2 m s _I over a 
5 min period before the Ethel, Texas, tornado of 
14 May 1974 and Shenk (1974) shows another example 
with a w of about 3 m s -1 . According to Shenk the 


mean diameter of the overshooting domes is approxi- 
mately 5.5 km which corresponds to a circular area 
of 24 km 2 . This area is an order of magnitude smaller 
than the appropriate satellite area of 15 2 or 225 km 2 . 
An example of w determined from radar top ascent 
can be determined from Burgess and Lemon (1976; 
see their Fig. 5.6). Prior to the Union City, Oklahoma, 
tornado, the supercell top ascended at ^3 m s -1 . This 
velocity is applicable to an area about 1-2 km on 
a side. 

From the preceding comparisons it is obvious that 
the vertical velocities at cloud top derived in this 
paper are smaller than those derived from other 
sources, but because they are applicable to an area 
much larger, they are still reasonable. Closer com- 
parison between the satellite information and radar, 
aircraft or other observations is necessary for full 
validation. 

8. Conclusions 

Digital infrared data from a geosynchronous satellite 
(SMS 2) have been used to study thunderstorm growth 
rates on 6 May 1975 in relation to the occurrence of 
tornadoes and other severe weather. Areal expansions 
of cold areas, delineated by blackbody temperature 
isotherms, are shown to be useful in monitoring 
thunderstorm growth rates, even when there is mini- 
mal decrease with time of the lowest cloud top 
temperature. 

An examination of five cloud elements having eight 
tornadoes clearly associated with them indicates that 
in seven of the eight cases the first report of the tornado 
took place during, or just after, a rapid expansion 
of cold areas indicating rapid ascent to the thunder- 
storm top on the scale observed by the satellite. The 
rate of ascent is estimated to be approximately 0.4 m s" 1 , 
averaged over an area of about 15 km on a side. 



VERTICAL VELOCITY, Wlms 1 ) 

Fig. 15. Mean vertical velocity and divergence for severe and 
non-severe thunderstorm elements. 
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Thirty-nine thunderstorm elements in the analysis 
area were compared in relation to various parameters 
derived from the satellite data and the occurrence of 
severe weather. Although the ground truth used for 
defining the occurrence of severe weather may have 
errors, in this study two satellite-based variables are 
correlated with tornadoes and hail. These variables 
are the minimum cloud top temperature during the 
lifetime of the element and the logarithmic rate of 
increase of cloud top cold area. 

The severe thunderstorm elements tend to be colder 
in terms of blackbody temperature and grow more 
rapidly than non-severe elements. For both variables 
the frequency distributions for severe and for non- 
severe elements were shown to be significantly dif- 
ferent. A discriminant analysis using both variables 
showed that on one side of the discrimination boundary 
four out of 23 cases had severe weather, while on the 
cold temperature, fast growth rate side 11 out of 
16 elements had reported severe weather. Therefore, 
based on this limited, dependent sample, these pa- 
rameters, which are indicative of intense convection, 
are correlated with severe weather on the ground. 
A potential warning lead time of 30 min was estimated 
by comparing times of severe weather reports and 
times the satellite information met the derived criteria. 

Equations were derived relating the areal expansion 
of blackbody temperature isotherms to vertical ve- 
locity and divergence. Mean vertical velocity profiles 
were calculated for severe and non-severe elements, 
with the results showing that the severe elements 
have a mean vertical velocity twice as large as the 
non-severe elements. The technique developed allows 
vertical velocities and divergence to be calculated 
from the satellite data on a horizontal spatial scale 
approaching that of individual thunderstorms. Infor- 
mation such as this will be useful in future work on 
tornado cloud top comparisons, further comparison 
of severe and non-severe elements, and for com- 
parison with and verification of various types of 
cumulus models. 

In summary, convection intensity information de- 
rived with short interval (^5 min) infrared digital 
data from a geosynchronous satellite can be statis- 
tically related to severe weather reports on the ground. 
Reasonable vertical velocity and divergence estimates 
can also be made from the data. The satellite data 
and the techniques developed have the potential to 
be a part of a severe thunderstorm detection and 
monitoring system. 
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ABSTRACT 

Quantitative observations of thunderstorms in the midwest United States made with short-interval 
(5 min) geosynchronous satellite data are examined in relation to concurrent digital radar observations 
for one case study over a limited area. Individual thunderstorms are defined in the satellite infrared (IR) 
data by the location of relative minima in the equivalent blackbody temperature ( T BB ) field. In a large 
majority of cases, these satellite-defined thunderstorms coincide with individual radar echoes. This agree- 
ment allows comparison of digital satellite and radar data for individual thunderstorms. 

The evolution of individual thunderstorms in terms of radar echo and satellite-observed cloud features 
is examined. An examination of a number of storms indicated that the first low-level radar echo (18 dBZ) 
appeared when the satellite observed cloud-top minimum T BB had a mean of 246 K (7.4 km). As the 
storms evolve, larger reflectivities appear as the cloud tops penetrate upward to colder temperatures. 

Larger reflectivity values (>50 dBZ) begin as the storms approach and penetrate the tropopause. 

Maximum radar reflectivity is shown to be correlated with satellite-based estimates of thunderstorm 
intensity. Thunderstorm top ascent rates in the 235-240 K (—8.8 km) region indicate the intensity of the 
initial storm updraft and are correlated with the maximum storm reflectivity with weak cells (~dT BB ldt 
of 1 K min' 1 ) having maximum reflectivity of 30-40 dBZ and strong cells (3-4 K min" 1 ) having echoes 
of 2*50 dBZ. The minimum T BB observed during the lifetime of the storm (r min ), indicative of maximum 
storm top height, is also correlated to maximum storm rainfall. Storms with tops colder (higher) than the 
tropopause (212 K) have the highest rainfall rates in the severe storm situation examined here. The 
parameter T min is also very well related to maximum volume rain rate as estimated from the radar data. 

Storms observed to reach temperatures lower than the tropopause temperature had volume rain rates 
of the order 10 3 m 3 s' 1 , compared to 10 2 m 3 s" 1 for weaker storms. 


1. Introduction 

In this study the evolution of thunderstorms is 
investigated using 5 min interval geosynchronous 
satellite data and concurrent digital radar data. The 
objective of the research is twofold: to determine the 
evolution of rainfall in individual thunderstorms in 
relation to satellite-observed cloud characteristics; 
and second, to determine the potential and limita- 
tions of using this type of satellite data to estimate 
convective precipitation on a storm-by-storm basis. 

The detection (and estimation) of convective rain- 
fall from visible and infrared (IR) image data is a 
complex problem; precipitation must be inferred 
from measurements of the radiative properties of the 
cloud top without any direct information on the pre- 
cipitation particles within the cells. Several schemes 
have been developed to use 30 min interval geosyn- 
chronous data to follow cloud masses and estimate 
their precipitation through empirical relations de- 
veloped through comparison of satellite-observed 
cloud parameters and radar-observed precipitation. 


1 A preliminary version of this paper was presented at the 

AMS Nineteenth Conference on Radar Meteorology, Miami, FL. 


Arkin (1979) has found correlations as high as 0.89 
between 6 h rainfall accumulation and the fractional 
coverage of the GATE B-scale array by clouds of 
varying (IR determined) heights. In more involved 
schemes Griffith et al. (1978) and Stout et al. (1979) 
have derived empirical relations between cloud area 
and echo area. Lovejoy and Austin (1979) found 
that a combination of IR and visible data was good 
for determining raining areas in convective regimes, 
but poor for determining rain rates, a conclusion also 
reached by Wylie (1979). 

The satellite-based precipitation estimation 
schemes already mentioned have primarily been 
tested over large areas (10 5 km 2 ) and appear to 
produce better results over longer time periods (e.g. , 
Woodley et al. y 1980). This latter technique follows 
the evolution of thunderstorms as depicted in the 30 
min data based on the outline of the 253 K isotherm. 
For the analysis of individual, rapidly changing 
thunderstorms higher time frequency is needed. 

In the current study, measurements of thunder- 
storm intensity determined from short-interval geo- 
synchronous digital satellite data are compared to 
precipitation measurements from digital radar 
data. Adler and Fenn (1979a) have shown that satel- 
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lite-based thunderstorm intensity estimates, such as 
the rate of blackbody temperature isotherm expan- 
sion and minimum cloud top temperature, are 
correlated with ground reports of severe weather. 
They also have found a significant difference in the 
mean calculated cloud top ascent rates between 
cloud elements associated with severe weather re- 
ports and those with no such reports (Adler and 
Fenn, 1979b). Severe weather reports represent iso- 
lated events, often not well located in space and 
time. The use of simultaneous digital radar informa- 
tion, remapped into the coordinate system of a 
geosynchronous satellite, allows a direct compari- 
son to be made between radar-observed precipita- 
tion and visible and infrared (IR) signatures at 
cloud top. Therefore, the evolution of precipitation 
in individual thunderstorms can be followed and 
compared to cloud evolution as observed with the 
satellite data. 

2. Data and methodology 

On 24 April 1975 intense, isolated convective 
storms developed in an area extending from central 
Oklahoma eastward to western Arkansas and 
Missouri during the late afternoon and early eve- 
ning. The synoptic and mesoscale analyses will not 
be presented here, but are available in Negri and 
Vonder Haar (1980) and Maddox et aL (1980). The 
radar data consisted of reflectivities at 0° elevation 
angle for WSR-57 radars at Monett, KiO; Kansas 
City, MO; and Oklahoma City, OK collected 
through the National Weather Service's Digitized 
Radar Experiment (D/ RADEX) program. The use of 
D/RADEX and its products is discussed in more 
detail by Saffle and Greene (1978) 2 . Areal resolution 
was 1.85 km by 2° of radar azimuth, and data were 
available nominally every 12 min in the form of nine 
discrete reflectivity levels. Conversion of radar 
reflectivity to rainfall rate was accomplished using 
the relation Z = 200/? 16 (Marshall and Palmer, 
1948), where Z is the radar reflectivity factor 
(mm 6 m~ 3 * * ) and R the rainfall rate (mm h _1 ). 

Utilizing the image processing capabilities of 
NASA’s Atmospheric and Oceanographic Informa- 
tion Processing System (AOIPS), D/RADEX reflec- 
tivities were transformed from polar (radar) coor- 
dinates to a Cartesian array, retaining a resolution 
of —10 km 2 and were subsequently remapped into 
the geosynchronous satellite’s coordinate system for 
direct comparison with cloud top features. The re- 
mapping was performed using bilinear interpolation 
with nearest neighbor resampling. The accuracy of 


2 Saffle, R. E., and D. R. Greene, 1978: The role of radar in the 

flash flood watch warning system: Johnstown examined. Pre- 

prints 1 8th Conf. Radar Meteorology, Atlanta, Amer. Meteor. 

Soc., 468-473. 


the remapped echo positions is -4 km and the 
areal extent of the remapped echoes is - 10% greater 
than the unremapped echoes, primarily in areas of 
low reflectivity. 

The satellite data set consisted of 5 min interval 
(rapid scan) imagery from 2200 GMT on 24 April 
to 0100 GMT on the 25th. Visible channel (0.6 fx m) 
radiances at 1 km resolution were used for naviga- 
tion and for the detection of overshooting domes 
against the cirrus anvil background. Quantitative 
use of the imagery involved the use of infrared 
(11.5 /xm) radiances at 8 km resolution. Radiances 
were converted to equivalent blackbody tempera- 
tures ( T bb ) and for each storm two parameters 
were noted: the rate of minimum cloud-top tempera- 
ture decrease (dT BB idt) and the minimum cloud- 
top temperature (T min ) observed during the thunder- 
storm’s lifetime. 

Examples of imagery used in this trispectral ap- 
proach to thunderstorm evolution are given in Figs. 

1 and 2. The first panel in each figure is a 1 km 
resolution visible image which delineates an area 
— 500 km on a side that encompassed most of the 
storms studied on this day. Fig. lb and 2b contain 
radar reflectivities from the three D/RADEX sites, 
combined and remapped into the perspective of the 
SMS/GOES satellite. The third and fourth panels are 
two different representations of the infrared images 
for these storms. Figs, lc and 2c have been en- 
hanced by the digital enhancement used opera- 
tionally by NOAA/NESS in the convective season. 
The definition of individual thunderstorms in the 
IR imagery often requires the contouring of every 
grey level (1 K) of temperature. Selected IR iso- 
therms based on this approach are presented in the 
final panel. 

The visible image at 2247 GMT (Fig. la) indicates 
three storms developing in south-central Oklahoma 
with several smaller storms forming on a flanking 
line to the west (refer to complex labeled C in 
Fig. la). As viewed in the IR (Figs, lc and Id) these 
main cells have corresponding T BB minima of 213, 
214 and 220 K, respectively. Subsequent calcula- 
tions are based on defining thunderstorm cells in 
terms of features often as small as those observed 
in Fig. lc (the smallest resolvable feature represents 
an area of —100 km 2 ). The corresponding radar 
display remapped into satellite coordinates (Fig. lb), 
indicates the precipitating portions of these storms, 
where grey represents reflectivities in excess of 
25 dBZ and black represents those in excess of 46 
dBZ or equivalent rainfall rates of 1 .3 and 27 mm h" 1 , 
respectively. In these young storms there is a one-to- 
one correspondence between the centers of maxi- 
mum precipitation and the T BB minima. An analysis 
of satellite IR features versus radar imagery for these 
storms also is given by Adler and Fenn (1979b). 
Other young storms (north of B) also show the 
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Fig. 1. (a) l km resolution visible image at 2247 GMT 24 April 1975. Area is —500 km on a side, (b) Reflectivity from radars 
at Monett, Missouri; Oklahoma City, Oklahoma; and Kansas City, Missouri; 2248-2252 GMT, remapped to match satellite 
geometry. Grey: 24-45 dBZ, black: >46 dBZ. (c) 8 km resolution 1R image enlarge to scale of visible image. Gray scale used is 
the digital enhancement used by NOAA-NESS (Mb curve), (d) Selected T BB isotherms from Fig. lc utilizing the one degree 
temperature resolution available. 


unique correspondence between individual echo and 
T bb minima. In a more mature storm (D) the rela- 
tion is less obvious with the coldest T BB (216 K) 
located downwind of both the intense portion of the 
radar echo and the cumuliform towers seen in the 
visible image. The area of 219 K temperature over 
the precipitation maximum may be an overestimate 


of cloud-top temperature due to the response of the 
IR sensor along the thermal gradient between the 
hot cloud-free environment and the cold cloud top. 

The coldest (207 K) and most intense storm at this 
time was the circularly shaped storm labeled A. The 
intense echo in this storm is located beneath a V- 
shaped pattern at cloud top (Fig. lc), possibly re- 
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Fig. 2. (a)-(d) Composite data as in Fig. 1, except 0012 GMT 25 April 1975. 


lated to blocking flow around the top of the storm. 
Other examples of this phenomena will be discussed 
subsequently. 

The later stages of these storms is illustrated in 
Fig. 2, the corresponding imagery 85 min later at 
0012 GMT on the 25th. The storms at C in Fig. 1 
have evolved into the line of convective towers X-Y 
(Figs. 2a and 2b), while those at B have expanded 


to form the line Y-Z. The storm at D has moved 
eastward and new convection has taken place to its 
north. Storm A has moved rapidly southeastward 
out of the area of interest. 

Between X and Y in Fig. 2, there is good agree- 
ment between the orientation of radar echoes and the 
line of overshooting tops in the visible image. How- 
ever, one of the two IR cold features (the larger 
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Fig. 3. Top: minimum blackbody temperature (K) versus time 
(GMT) for cloud 24. Bottom: Echo area (km 2 ) at selected dBZ 
thresholds versus time (GMT) for cloud 24. 


one) at 206 K is centered north of the X-Y axis, 
thereby giving a false impression of storm alignment. 
The storm associated with the large 206 K area has 
passed its maximum intensity, and the center of the 
cold area may not be colocated with the remaining 
updraft and rainfall. To the southeast of this area, 
a large, enclosed area of 209 and 210 K is also not 
apparently associated with any precipitation. The 
visible channel is particularly useful here in delineat- 
ing these minima as smooth anvil material, however, 
it is evident that for mature storms the location of 
intense precipitation and the apportionment of pre- 
cipitation by satellite estimation schemes cannot be 
based solely on infrared observations of the cirrus 
anvil. 

Often, storms cannot be observed by the satellite 
until they penetrate the cirrus anvil generated by 
previous convection. Such is the case with the line 
of tops Y-Z, which is spatially correlated with the 


line of echoes and defined by two T BB minima of 
206 and 210 K, both colder than the tropopause 
temperature and in this case are also colocated with 
the twin intense echo maxima. The 206 K storm be- 
tween Y and Z also exhibits a V-notch pattern in the 
T bb field frequently associated with intense thunder- 
storms (Adler et al ., 1981). This type of feature is 
probably related to obstacle flow at anvil level 
around a strong, erect updraft. There is much evi- 
dence for such an effect (Browning and Foote, 
1976). The warm (in terms of T BB ) area downwind 
is probably caused by the subsiding flow to the rear 
of the storm top, although there is a possibility it is 
related to the presence of lower stratospheric cirrus 
clouds formed downwind of the overshooting top. 
If these clouds were present, they would appear 
warmer than the anvil background. 

In this study individual thunderstorms were de- 
fined in the satellite data using a sequence of IR 
and visible images and a technique described by 
Adler and Fenn ( 1979a,b). Points of relative minimum 
in T bb were tracked and defined as thunderstorms 
if they displayed good temporal continuity, a cumuli- 
form appearance in the visible data, and a period of 
rapid change in the cloud top T BB . These subjective 
parameters are indicators of deep convection and 
eliminate anvil debris and other anvil perturbations 
from consideration as thunderstorms. The thunder- 
storms so defined in the satellite data agreed very 
well with individual radar echoes, especially in the 
growing and maximum rainfall rate stages. During 
later stages, when the light rain areas on the 
periphery of the echoes merge as the storms form 
into lines, the identified IR features still are well 
associated with radar reflectivity maxima, although 
they are often not as well colocated as during 
earlier stages. Infrequent (in this study) mergers of 
storms are handled by defining the merged storms 
to be a continuation of most intense of the two 
previous storms. The defined lifetime of the weaker 
storm ends at the time of merger. 

The high time-resolution satellite data are im- 
portant in the identification and tracking of individual 
thunderstorms. From a solitary image, or from 
coarse time resolution (e.g., 30 min) data, it would 
become very difficult to eliminate T BB anvil features 
hot associated with thunderstorm updrafts and it 
would also be more difficult to determine whether 
features on successive images, in fact, are manifes- 
tations of the same thunderstorm. 

3. Evolution of individual thunderstorms 

In this section three examples of evolving thunder- 
storms will be discussed with respect to satellite 
infrared and digital radar observations. The results 
of the analysis of a weak storm are shown in Fig. 3. 
The specific thunderstorm, cloud 24, is evident in 
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Fig. 1 in the extreme southwest corner of the 
images with a cloud top T BB of 240 K at 2247 GMT. 
The growing thunderstorm is noted in the IR data as 
an ascending cloud top (decreasing IR Tbb values) 
going from 260 to 222 K over a 40 min period (Fig. 
3). The time rate of T BB decrease (-dT BB /dt) was 
approximately 1 K min -1 , a small value indicating a 
weak storm (Adler and Fenn, 1979b). After 2315 
GMT the storm could no longer be identified in the 
IR data due to anvil cirrus from a stronger storm. 
The first echo (18 dBZ) appeared at 2236 GMT at 
a satellite-observed cloud-top temperature of 254 K 
(Fig. 3). Echo area expanded sharply and the maxi- 
mum reflectivity increased from 2230 to 2300 GMT. 
Then a slower decrease began until the echo essen- 
tially disappeared at approximately 0000 GMT. 
Cloud 24 had a small, relatively weak echo and 
also appeared weak in the sequence of satellite 
digital data with slow (1 K min -1 ) cloud top ascent 
rates and no observed temperatures lower than 
222 K, 10° above tropopause temperature. Some- 
what lower temperatures may possibly have been 
observed if the cloud top had not become blocked 
from the view of the satellite by anvil cirrus from 
an adjacent, stronger thunderstorm, but the mini- 
mum observed T B b is coincident with the radar- 
observed maximum intensity. 

A second example is cloud 22, an intense hail- 
storm (7 cm diameter hail), which is evident in Fig. 
1 as the easternmost cell in area C with a T BB of 
213 K at 2247 GMT. Cloud 22 was a very intense 
thunderstorm with values of -dT BB !dt of 3.5 K min -1 
in the 235-240 K region (see Fig. 4a), putting it into 
the severe thunderstorm range (Adler and Fenn, 
1979b). The first radar echo (18 dBZ) appeared at a 
T bb of 238 K at 2236 GMT, increased rapidly in 
magnitude to a reflectivity of 46 dBZ at 2248 GMT 
and to a storm maximum of 50 dBZ at 2300 GMT 
(Fig. 4b). Maximum reflectivity decreased after 2315 
GMT. The satellite data indicate a thunderstorm 
top rapidly ascending through the upper troposphere 
until the tropopause (212 K) is approached and pene- 
trated. Expansion of T BB isotherms in the anvil is 
displayed in Fig. 4c. During the period of rapid echo 
expansion and intensification (2240-2300 GMT) 
there is also a decrease in minimum T BB and an 
expansion of T BB isotherms in the storm’s anvil. 
Nearly simultaneous with the maximum intensity of 
the radar echo (2300-2312 GMT) is the satellite 
observed highest (coldest) cloud top (206 K) at 
2316 GMT. After this time the echo magnitude 
drops and by 2335 GMT the distinct T BB relative 
minimum associated with cloud 22 has disappeared 
into the broad background of the anvil. This would 
seem to indicate a weakening of the updraft. 
Significant rain, although not the heaviest, continues 
to fall even after the storm top sinks back into the 
cirrus background. This may be a persistent problem 


in any attempt to relate storm-scale satellite ob- 
servations with convective rainfall in the mature 
phase of storms. On the other hand, the heavy rain 
did not begin until about the time the tropopause 
was reached by the cloud. This may be important 
because many storms cannot be observed until they 
penetrate the cirrus shield produced by previous 
convection (Adler and Fenn, 1979a). Thus, even 
these storms, only observable in their mature stage, 
can be identified near the time the heavy rain begins. 

A third example, cloud 19, was a very intense 
storm in its mature stage, as can be seen in Fig. 2, 
where it is located in extreme northwest Arkansas. 
However, as shown in Fig. 5a, it began as a weak 
storm with a -dJ BB ldt of approximately 1 K min -1 
in the 235-240 K height region. During the early 
portion of the storm, the radar echo also remained 
weak in agreement with the slow cloud-top ascent 
rate. At approximately 2300 GMT, the character of 
the thunderstorm as viewed by both the satellite 
and the radar changed, with a rapid intensification 
of the echo and a rapid expansion of the anvil 
edge and the T BB contours within it (Figs. 5b and 
5c, respectively). The reason for this abrupt change 
from weak to intense storm is unclear, but may be 
due to a merger of updrafts, an enhancement of 
the low-level convergence through storm superposi- 
tion over a low-level boundary, or other effect. The 
cloud temperature continues to decrease during this 
rapid echo growth period and finally reaches a 
minimum temperature of 204 K (~8 K below tropo- 
pause temperature), slightly before the appearance 
of the storm’s maximum reflectivity, 62 dBZ. 

In summary, the three example storms show the 
range of thunderstorms analyzed in this case study. 
Cloud 24 is a weak, short-lived storm with a small 
satellite-observed growth rate and a resultant small 
maximum height (minimum temperature). Moderate 
rain is only apparent for a short time and the entire 
life cycle of the radar echo is <90 min. The small 
size of the echo and the moderate rain values com- 
bine to produce a small volume rain rate. Cloud 
22 is an intense, rapidly evolving severe thunder- 
storm with a large cloud-top ascent rate in the 
upper troposphere, a low minimum T BBi and a large 
maximum reflectivity. The period of heaviest rain 
(50 dBZ) was relatively short (20 min). Cloud 19 
evolved into one of the heaviest raining storms on 
this day, despite appearing initially as a slow grow- 
ing, weak storm in both the radar and satellite 
data. The heaviest rain in cloud 19 was long-lived 
producing an enormous total storm rainfall. 

4. Results for all storms 

A total of 15 thunderstorms were analyzed with 
both the radar and satellite data on this day. While 
a majority of storms were observed in the data 
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throughout their entire lifetime, a few could only be 
identified after they penetrated a cirrus cloud deck. 
Other weaker cells became obscured during their 
later stages by anvil cirrus from other storms. 

Fig. 6 is a scatter diagram showing the thunder- 
storm top T B b value at the time when the storm’s 
reflectivity initially reached a particular magnitude. 
The correlation coefficient between these two pa- 
rameters was -0.85. First echo (18 dBZ) was noted 
between 229 and 260 K with a mean temperature 
of 246 K. This mean temperature is equivalent to a 
height of -7.4 km. Browning and Atlas (1965) and 
Atlas (1966) have examined this question of first 
echo height observationally and theoretically, 
respectively. The former study found the midpoint 
of the first echo in a Oklahoma severe storm to be 
at or above the -30°C (243 K) level, significantly 
colder than the similar level (-0°C) in non-severe 
convective clouds (Battan, 1963). Although the re- 
gression line in Fig. 6 also points to a first-echo 
height relationship, the results are not strictly com- 
parable. As previously mentioned, the IR measure- 
ments tend to overestimate the true cloud tempera- 
ture due to field of view and response problems. 
In addition, the D/RADEX data were available only 
at 0° elevation and 12 min intervals and archived 
with a minimum reflectivity of 18 dBZ. Thus, the 
derivation of a first-echo height relationship in the 
manner of Browning and Atlas (1965) was not 
possible. 

As the storms evolve, greater reflectivities 
appear as the storms penetrate the upper troposphere 
to colder temperatures. Larger reflectivity values 
(>46 dBZ) begin as the storms approach and pene- 
trate the tropopause. The large scatter of points 
along the T BB axis at the warm temperature end of 
Fig. 6 is partly due to the time resolution of the 
radar data (12 min). With rates of change of T BB of 
2 K min -1 , a time error of 6 min in the appearance 
of a reflectivity value will translate into a 12 K error. 
At colder temperatures, where ~dJ BB ldt is smaller, 
the radar time resolution is less critical and the 
magnitude of the scatter decreases. 

Adler and Fenn (1979b) have shown that one 
measure of thunderstorm intensity is the cloud top 
ascent rate (- dT BB /dt ) of storms in the 235-240 K 
(8.8 km) region. Storms with large ascent rates 
(>2 K min -1 ) tend to have associated severe weather 
reports. Fig. 7 shows the relation of this satellite 
based intensity value and the maximum D/RADEX 
reflectivity, where each data point is an individual 
(numbered) thunderstorm. The reflectivity value 
used in the diagram must have existed in the storm 
for three successive scans (>24 min); this eliminates 
short-lived higher values. For those storms observed 
by the satellite in their immature stage, there is a 
weak correlation (-0.66) between cloud-top ascent 
rate (proportional to updraft magnitude) and maxi- 


mum reflectivity (Fig. 7), however, this correlation 
becomes 0.9 if the outlier (cloud 19) is excluded. 
This was an intense thunderstorm which began as an 
apparently weak cell and whose satellite and radar 
growth curves were presented in the preceding 
section. 

The four storms viewed by the Oklahoma City 
radar with growth rates of <2 K min -1 are clearly 
weak thunderstorms (clouds 24, 30, 31 and 32). These 
storms also had relatively light rain (as indicated by 
low maximum reflectivity). The more intense storms 
had progressively larger maximum reflectivity as the 
satellite parameter varied from 2 to 4 K min -1 . 

Another satellite parameter that indicates storm 
severity (Adler and Fenn, 1979a) is the minimum 
blackbody temperature (T min ) achieved by the storm 
during its lifetime, observed by following an in- 
dividual thunderstorm with rapid interval IR imagery. 
A plot of J min versus the corresponding observed 
maximum radar reflectivity for the 15 identifiable 
storms is presented in Fig. 8. Again the reflectivity 
maxima had to exist for at least 24 min. In addi- 
tion, the occurrence of J min and maximum re- 
flectivity was not necessarily simultaneous, but oc- 
curred on the average within 15 min of one another. 
Despite the small sample size, there is an excellent 
correlation (-0.92) between T min and maximum re- 
flectivity. In addition, all storms observed to pene- 
trate the tropopause (212 K) had reflectivities in 
excess of 50 dBZ. However, the relationship of 
rainfall intensity to storm intensity on the heavy rain 
side is probably more complicated and is influenced 
by such things as precipitation efficiency, entrain- 
ment, and interactions with the environment and 
neighboring storms. 

A parameter particularly useful in identifying the 
flash-flood-producing thunderstorm is the volu- 
metric rain rate (VRR). The VRR is computed from 
the radar data as follows: For every isolated radar 
echo, the Marshall-Palmer relation Z = 200/? 16 is 
applied to each reflectivity value to obtain the 
equivalent rain rate (mm h -1 ). The product of that 
rain rate with its areal extent (km 2 ) is then formed 
and summed over the nine possible reflectivity levels 
present in a storm. When echoes merge, the pro- 
cedure remains the same, except that a subjective 
determination of the boundaries of the defined rain 
area is made with the point of maximum reflectivity 
remaining well inside the defined boundaries. This 
does not introduce significant bias because the lower 
reflectivities that exist on the edges of these sub- 
jectively defined echoes have small equivalent rain 
rates (0.5-10 mm h -1 ) and hence comprise only a 
small percentage of the total VRR. 

The comparison between VRR (m 3 s -1 ), plotted 
on a logarithmic scale, and minimum T BB (Fig. 9) 
reveals a remarkably linear relationship. Each data 
point represents an individual thunderstorm defined 
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by a subjectively defined area centered on a reflec- 
tivity maximum and a T titi minimum determined to 
be associated with active convection. In general, 
the maximum VRR occurred nearly simultaneously 
with the maximum reflectivity. Storms that were ob- 
served to penetrate the tropopause had VRR values 


that, in general, were an order of magnitude larger 
than the warmer (lower) storms. The extremely good 
correlation (-0.95) suggests that we are relating two 
meaningful parameters in a physically sensible way. 
The T min reached in each storm is a measure of the 
maximum updraft intensity averaged over a substan- 
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Fig. 6. Satellite-observed, thunderstorm top T BB at the time of first appearance of various 
radar reflectivity values for all storms observed on 24-25 April 1975. 



Fig. 7. Maximum observed reflectivity of 24 min duration 
(dBZ) versus cloud-top ascent rate (K min" 1 ). Each data point 
is an individual thunderstorm defined in both the IR and radar 
data. Scale on right shows equivalent rainfall rate (mm h" 1 ) using 
the relation Z = 200 R 1 6 . 


tial area 100 km 2 ), while the VRR is an integrated 
measure of the rainfall, because rain generated in 
the updraft may extend into areas beyond the updraft 
itself. These results are particularly interesting be- 
cause they not only provide evidence of a very plau- 
sible physical relationship (to be addressed in subse- 
quent research) but, if validated by further observa- 
tions, provide a basis for eventual thunderstorm 
rainfall estimation from satellite data using a physi- 
cally based method. 

The flash flood potential of a thunderstorm through 
accumulation involves not only its instantaneous 
rainfall rate, but also its rate of motion. A compari- 
son of the radar determined, accumulated rainfall 
for two storms with similar reflectivity patterns was 
undertaken. These storms are cloud 4, the large cir- 
cularly shaped storm north of “D” in Fig. 1 and 
cloud 19, previously discussed, located in northwest 
Arkansas in Fig. 2. Cloud 4, whose echo centroid 
traveled at 18.8 m s _I , left a narrow elongated swath 
of precipitation 10-20 mm deep across central Mis- 
souri, reaching a maximum of 28 mm. Cloud 19, how- 
ever, moved somewhat slower (12.8 m s -1 )> which 
allowed accumulations of 65 mm of rain along the 
Oklahoma-Arkansas border during the same time 
span. A similar difference in storm speeds is ob- 
tained by following the movement of the correspond- 
ing T bb minima in the satellite imagery. Thus the 
use of the IR data, applied on a thunderstorm scale, 
allows for not only the inference of storm intensity 
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Fig. 8. Maximum observed reflectivity of 24 min duration (dBZ) versus minimum blackbody 
temperature. Note high reflectivities associated with tropopause penetrating storms. 


but also the rate of storm motion. In addition, a 
new thunderstorm (cloud 31a) developed to the west 
of cloud 19 and by 0124 GMT maximum radar de- 
termined accumulated precipitation in parts of west- 
ern Arkansas had reached 72 mm. 

5. Conclusions 

The purpose of this study was to investigate the 
evolution of thunderstorms of various intensities 
with short interval (5 min) geosynchronous digital 
satellite data and with digital radar data (12 min 
interval) to determine the potential and limitations 
of using geosynchronous satellite data to detect re- 
gions of convective precipitation and estimate the 
rainfall intensity. Thunderstorms on 24 April 1975, 
in eastern Oklahoma and surrounding areas were 
examined with SMS/GOES digital satellite data and 
D/RADEX digital radar data from the Oklahoma 
City, Oklahoma, Monett, Missouri, and Kansas City, 
Missouri radar sites. The radar data were remapped 
into satellite coordinates to facilitate display and 
comparison with the satellite data. 


Individual thunderstorms were defined in the satel- 
lite IR data by location of relative minima in the 
equivalent blackbody temperature (T BB ) field. In 
virtually all cases, the locations of these satellite- 
defined thunderstorms coincided with individual 
radar echoes, notably in the early stages of develop- 
ment. This agreement allows comparison of digital 
satellite and radar data for individual thunderstorms . 

Maximum radar reflectivity was shown to be cor- 
related with satellite based estimates of thunderstorm 
intensity. Thunderstorm top ascent rates in the 235- 
240 K (8.8 km) region reflect intensity of the storm 
updraft and are correlated with the maximum storm 
reflectivity with weak cells ( -dT BB fdt of 1 K min -1 ) 
having maximum reflectivity of 30-40 dBZ (3-5 mm 
h -1 ) and strong cells (3-4 K min -1 ) having echoes 
of >50 dBZ (49 mm h -1 ). 

For mature storms, the maximum updraft intensity 
is estimated (in a relative sense) by the minimum 
T bb observed during the lifetime of the storm (r min ). 
This satellite-based parameter is inversely related to 
maximum cloud height and was shown in this study 
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to be correlated to maximum storm radar reflectivity. 
Storms with 7 mln colder (higher) than the tropopause 
had the highest radar reflectivities. The parameter 
T min also is very well related to maximum volume 
rain rate as estimated from the radar data. Storms 
observed to reach temperatures lower than the trop- 
opause temperature had volume rates of the order 
10 3 m 3 s _1 . A physical basis for these statistical 
relationships between the satellite-observed param- 
eters and radar-estimated rainfall is currently being 
investigated. 

The first low-level radar echo (18 dBZ) appeared 
when the satellite-observed cloud-top minimum T bb 


had a mean of 246 K (7.4 km). As the storms evolve, 
larger reflectivities appear as the cloud tops pene- 
trate upward to colder temperatures. Larger reflec- 
tivity values (>50 dBZ) began as the storms ap- 
proached and penetrated the tropopause (212 K). 

These results indicate that inferences on thunder- 
storm updraft intensity can be made from such 
satellite derived quantities as dT HB idt and T min . For 
individual thunderstorms on 24 April 1975, such in- 
tensity measurements were related to storm precipi- 
tation rate as indicated by concurrent digital radar 
data. The conclusions support the effort to use satel- 
lite data to detect heavy convective precipitation 
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and estimate its magnitude, but emphasize the need 
for high time-resolution (3-5 min) data. The short 
interval data are necessary in order to 1) accurately 
define and follow individual thunderstorms unam- 
biguously; 2) accurately measure rate of change pa- 
rameters on a time scale appropriate for thunder- 
storms; and 3) accurately determine the minimum 
T bb achieved during the thunderstorm’s lifetime. It 
must also be emphasized that the empirical relations 
between T BB parameters and rainfall rate derived 
here do not apply to all data points on an IR image, 
but only to the locations of defined thunderstorms 
in that image. In addition, we have examined only 
one day and only one area. These relationships may 
be qualitatively similar but quantitatively different 
elsewhere, notably in tropical regions. 
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I . INTRODUCTION 

Estimation of convective rainfall from 
visible and window infrared (IR) satellite data 
continues to be of great interest, especially for 
heavy convective rain. While radar techniques 
provide an estimate of rainfall rate based on 
physical relations between the radar return 
signal and the cloud rainfall itself, the existing 
satellite-based IR and visible methods (e.g., 
Griffith et al., 1978, 1981; Scofield and Oliver, 
1977; Stout et al., 1979) require rainfall rate 
information to be inferred from cloud top features 
and equivalent blackbody temperature (T RR ) data. 

The cloud/rainfall relationships used are deter- 
mined empirically. The relations are derived for, 
and in turn applied to, time and space scales much 
larger than that associated with individual thunder- 
storms. One of the major drawbacks to such tech- 
niques is the lack of a strong physical basis, 
which makes it difficult to account for geographic, 
seasonal, and other variations. 

Recently, there have been attempts to 
estimate the moisture budget (Park and Sikdar, 

1980) and rainfall (Mack and Wylie, 1981) in 
large convective complexes by estimating anvil 
mass outflow from satellite data in order to cal- 
culate convective vertical mass transports. This 
information, in addition to moisture data (from 
rawinsonde) and an assumed entrainment factor, 
leads to an estimate of precipitation in the 
complex based on a simple model of cumulus mass 
and moisture transport discussed by Austin and 
Houze (1973). This type of approach has the 
potential to make physically-based estimates of 
large scale rainfal 1 . 

Major questions remain as to what is 
occurring on the scale of an individual thunder- 
storm and whether or not we can estimate any pre- 
cipitation rate parameters on that scale. Negri 
and Adler (1981) have examined an admittedly small 
sample (15) of thunderstorms in Oklahoma and 
surroundinq states on 1 day with short interval 
(5 minutes) GOES data and digital radar data. 
Individual thunderstorms were defined in the satel- 
lite IK data by the locations of relative minima in 
the equivalent blackbody temperature (T PR ) field. 

The feature must also maintain continuity, appear 


convective In the visible data and undergo a 
period of rapid change in T RR In order to be 
defined as a thunderstorm, in virtually all 
cases, the locations of these satellite-defined 
thunderstorms coincided with individual radar 
echoes, notably in the early stages of develop- 
ment. This agreement allows comparison of digital 
satellite and radar data for individual thunder- 
storms . 


Maximum storm rainfall rate and maxi- 
mum storm volume rainfall rate (VRR) were shown 
to be correlated with two satellite-based esti- 


mates of relative updraft intensity: the cloud 

top ascent rate (-dT BR /dt) in the 235-240K 
(8.8 km) region, and°the minimum T RR (T . ) 
achieved during the storm’s lifetime. m1n 


The objective of this paper is to 
derive some simple expressions relating rainfall 
rate to thunderstorm updraft magnitude which can 
be interpreted in terms of the satellite observ- 
ables. These expressions will then be compared 
to the observations presented by Negri and Adler 
(1981), hereafter referred to as NA. The hope is 
to establish the gross basis for physical rela- 
tionships underlying the statistical results. 
Comparison of the satellite and radar observations 
with one-dimensional model results is also made. 


2. RAINFALL RATE EQUATIONS 

The process of rain formation in a 
thunderstorm is a complex interaction of storm 
dynamics and microphysics. However, the empirical 
results of NA indicate that at least for the one 
case study, the variation in maximum rainfall 
rate and maximum volume rainfall rate among 
storms is strongly related to the relative 
updraft intensity as indicated by the satellite- 
observed cloud top ascent rate or by T . . The 
attempt in the following paragraphs wili n be to 
derive pertinent equations which can be evaluated 
with the satellite observations. 


The total condensation rate C can be 
related to updraft intensity and moisture by start- 
ing with the moisture conservation equation. 


121 


assuming steady state conditions and horizontal 
uniformity and vertically Integrating the equa- 
tion from cloud base to cloud top. 


C 



PW dzda. 


(1) 


where p is the density, w is vertical velocity, 
q is saturated specific humidity, and z* and z fe 
a?e cloud top and base, respectively. Tne area 
A is the updraft area. Eq. (1) can also be 
expressed as 

C = - { I pwdadq, (2) 


where q is the saturated q at cloud base (also 
equal t8 surface q, assuming unmixed ascent to 
cloud base) and zero is the approximate cloud top 
saturated q. Of course, cloud top q can be greater 
than zero, which should be kept in mind for gen- 
erality. But for the data set of NA, the weakest, 
or shallowest, storm had a T . of 231K, or a 
cloud top q of 0.4 g/kg, conpired to a q of 
15 g/kg. Even for a cloud with a maximum top of 
only about 7 km, the cloud top q would be approxi- 
mately 10% of the cloud base value. Therefore, the 
assumption of q=0 at cloud top is acceptable for 
this exercise, but would have to be modified for 
shallow convection. 


Since we are dealing with saturated 
processes in the cloud, all q's will be saturated 
and the subscript will be dropped for convenience. 
The density p is assumed to be horizontally uniform 
and can be removed from inside the area integral. 
The updraft area integral 

| wda = wA (3) 

A 


where w is the average w over the updraft area. 

The satellite observations, through the cloud top 
ascent rate, give us an estimate of an area- 
averaged w, but not simply averaged over the up- 
draft area, but over the satellite IR Instantaneous 
Field of View (IFOV), approximately 10 km by 10 km. 
However, if w=0 outside the updraft and the satel- 
lite IFOV is always larger than the updraft area, 
then 


j wda = wA = w A , (A <A ) (4) 

A s s s 


where w e is the average vertical 

- . s - 


area (A*) of the satellite IFOV. 
becomes 5 


C = -A. 


pw s dq 


velocity over the 
Therefore, Eq. (2) 


(5) 


The vertical coordinate q (actually 
q ) is monotonic and has the characteristic that 
p s is approximately linear in q so that p can be 
expressed as: 

p = a o + a l 


where a and a, are constants appropriate for a 
particular mom adiabat e . The linear approxi- 
mation is valid because both variables decrease 
exponentially with height. Figure 1 shows the 
relationship of p, q and height along a moist 
adiabat (e =348K) ana the linear approximation 
used. Along a moist adiabat p is only a function 
of q and a more complete formulation could use 
the Derived relation. 



Figure J. Vertical profile of air density p as 
a function of a along the $ e =348K moist adiabat 
(solid line) ar& the linear approximation 
(dashed line) used in the calculations. 


The variable w (the subscript s is 
dropped for convenience) is parameterized in 
terms of maximum vertical velocity, w max , so that 

w = W max (1 - B q/ V (7) 

This formulation gives w=w at cloud top and 
a linear profile (with respect to q ) down to 
w=w (1-3) at cloud base, where 0<B<1 . The 
parSftiter B will be determined by an additional 
constraint. A linear profile of w with respect 
to q is a reasonable approximation as can be seen, 
for ixample, in the results of a one-dimensional 
cloud model run (Fig. 2). The w profile given by 
Eq. (7) does not include the decrease of w with 
height in the outflow layer at the top of the 
storm. However, this height region contributes 
very little to the condensation rate because of 
its low moisture values. 

Substituting Eqs. (6) and (7) into 
Eq. (.5) and integrating produces: 

C - [0 - e/2)a 0 + (0.5 - P/3)a 1 ]w max q 0 A $ (8) 

In this simple formulation, the total condensation 
rate, C, is proportional to the product of maximum 
updraft intensity and cloud base saturation spefi- 
fic humidity (also equal to surface specific 
humidity) . 
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Height above Cloud Base (km) 





Figure 2. Vertical profile of one- dimensional 
model vertical velocity plotted as a function 
of q along the Q & -348K adiabat . The model 
velocity is nearly linear with respect to q . 


In order to determine 6, we will assume 
that the only moisture source for the cloud is 
through cloud base. Since the system is assumed 
to be steady state, the flux of moisture through 
cloud base must equal the condensation rate, i.e.. 


C = Wo A s 

(9) 

where the subscript denotes cloud base 
From Eqs. (6) and (7) 

values . 

” ci + a, and 
0 0 l 

(10) 

w o * w nax (1 ” 8) 

(11) 

and Eq. (9) becomes 


C = (a o + “l )w max (1 ' e)q o A s 

(12) 


Solving tqs. (8) and (12) for 3 gives 


* = cx Q + a’ 4/3 (13) 

For the April 24, 1975, case study 
appropriate constants for the density profile 
are a = 35xl0 3 grrr 3 and a,=.69xl0 3 gm . Sub- 
stituting these values into Eq. (10) results in 
3= .54, which means that cloud base vertical 
velocity w =.46 w . At first, a vertical 
velocity profile with the cloud base value al- 
ready nearly one-half of the maximum value seems 
unreasonable and would be if we were considering 
peak velocity values at each level. However, one 
must remember that we are dealing with an average 
value over a horizontal area larger than the up- 
draft (the satellite IFOV ^100 km 2 ). In terms 
of vertical mass flux, this vertical velocity 
profile will give approximately constant mass 
flux with height, reasonable in our no entrainment 
model. For example, a w of 4 ms -1 (an inter- 
mediate value) corresponos x to a w of 1.8 ms" . 
Using an area of 100 km 2 and p ot 10 3 gm’ 3 re- 
sults in a cloud base mass flux estimate of 
1.8X10 11 g s’ 1 , similar to an average value for 
hailstorms (2.3X10 11 g s’ 1 ) calculated by Auer 
and Marwitz (1968) from aircraft data. 


Substituting the values for a , a, and 
6 into Eq. (8) gives 

C = .475 x 10 3 Wmax q 0 A $ (14) 

where the constant has units of gm . 

The volume rainfall rate (VRR) will be 
less than the total condensation rate, C, due to 
various loss terms. This is expressed simply as 

VRR = — n 5) 

p w 

where p is the density of water and E is the 
storm precipitation efficiency (0<^E<_1 ) . Therefore, 

VRR = .475x10” 3 E A $ q Q (16) 

Eq. (16) indicates that VRR is propor- 
tional to w . Since the area (A ) is constant, 
this actual rymea ns that VRR is proportional to 
vertical mass flux. Cloud top ascent rates in 
the 235-240K (8.8 km) region as used by NA can be 
converted to estimates of equivalent w^ through 
an appropriate lapse rate (Adler and Fenn, 1979). 

In this case, a lapse rate of 8.3 km r 1 was used 
which converts a cloud top ascent rate of IK min’ 1 
to a w TOx of 2 ms’ 1 . 

The 235-240K layer is used because it 
is intermediate to the warmer layers where the 
clouds are difficult to define and the colder 
layers where the ascent rates are affected by 
the tropopause. It is the layer with generally 
the largest ascent rates (Adler and Fenn, 1979). 

When we use the satellite cloud top ascent rates 
as a substitute for w in Eq. (16) we are 
assuming that the ascent rate during the cloud 
growth stage is equal to the maximum vertical 
velocity (area-averaged) In the interior of the 

storm at a later time. This is -probably not 
true, but the two parameters are correlated and 
the cloud ascent rates produce reasonable vertical 
mass fluxes (see previous paragraph). 

Figure 3 shows the results of calculations 
with Eq. (16) with q =15 g/kg (a reasonable value 
for the case study) and E=1.0. The data points 
represent the eight storms in the case study (NA) 
where the cloud top ascent rate in the 235-240K 
layer could be observed along with the maximum 
radar-observed VRR (instantaneous) during the 
storm's lifetime. The calculated line follows the 
data reasonably well, although the more intense 
storms straddle the line while the weak storms fall 
significantly below the line. The instantaneous 
maximum VRR may overemphasize transient peaks in 
the rainfall. An average over a reasonable time 
period may make a better comparison with our cal- 
culated values. Therefore, a 30 minute average 
(3 radar observations) was calculated and plotted 
in Fig. 4. The calculated VRR's have been reduced 
by 15-20% by the time averaging. 

The plots in Figs. 3 and 4, especially 
Fig. 4, indicate that the linear relation in 
Eq. (16) is a reasonable first approximation and 
that, for these storms, the VRR is roughly pro- 
portional to the vertical mass flux. Data points 
below the calculated line imply efficiencies (E) 
less than one. For the four intense storms, the 
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E value calculated from Fig. 4 is 80%, while the 
weak storms average 25%. Although the absolute 
value of the efficiencies should be treated very 
carefully because of all the assumptions in the 
derivation of the theoretical VRR and in the radar 
and satellite data analysis, the trend toward 
higher efficiencies for more Intense, larger 
storms is realistic because of the relation of 
entrainment to storm or plume size (Turner, 1962). 



W M *x lm« 'l 


Figure 3. Volume rainfall rate (VRR) calculated 
as a function of area-average from Eq (16). 

Dots represent observations fromNA of maximum 
instantaneous VRR and satellite-observed cloud 
top ascent rates. 



Figure 4. Same as Fig. 3 except VRR observations 
(dots) are time averages over SO minutes. 


Often a storm's cloud top ascent rate in 
the 9 km region cannot be measured because of anvil 
cirrus produced by previous convection. However, 

T - , the minimum temperature achieved during a 
s?i!Vs lifetime, can usually be observed, 
especially in intense storms. Because this para- 
meter is a relative measure of maximum storm 


height, it is also a proxy variable for w . A - 
relation between T . and w can possibly be 
derived through analysis of S one-dimensional cloud 
model (see section 3), or empirically derived as 
depicted in Fig. 5, where a curve is fitted to the 
satellite-based w and T , data for eight 
storms. The empirical relation 


W 5 - 32x10 ' 3(237 - T n,1n )2+ - 08 


(17) 


can then be substituted into Eq. (16) to obtain 

VRR=4.75x10‘ 4 [5.32x10‘ 3 (237-T . ) 2 +.08]E q A $ (18) 


The curve is plotted in Fig. (6) along with the 30 
minute average VRR data. Again, the weak storms fall 
considerably below the curve while the intense storms 
scatter about the curve. Points falling above 
the curve imply efficiencies greater than 1.0, 
which are only possible if there is an additional 
source of moisture above cloud base and there are 
no moisture losses, such as outflow in the anvil. 



Figure 5. Empirical relation of satellite- 

observed cloud top ascent rates to satellite- 

observed minimum T nrt > ., for storms in NA. 
du mzn 


The calculated curve, Eq. (18), gen- 
erally follows the trend of the data and the 
increase in slope as lower T . values are 
approached. However, it doesnot appear to cap- 
ture the sharp change in slope around the tropo- 
pause temperature (214K). This is partially 
because the empirical relation between w and 
T min (Eq. (17)) smooths out the transition. 

The apparent trend in efficiency is 
again present in Fig. 6. One term affecting 
efficiency is entrainment, and entrainment is 
usually parameterized as inversely proportional 
to cloud horizontal size. Since cloud width is 
related to cloud depth (as indicated by I n i n ), 
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we should be able to parameterize the entrainment 
and the efficiency in terms of cloud size. A modi- 
fication to Eqs. (8) and (18) incorporating an 
entrainment parameterization has been derived, but 
is not presented here due to lack of space. The 
modification reduces the calculated rainfall for 
the small clouds and tends to bring the calculated 
values into better relative agreement with the 
observati ons. 



Figure 6 . Volume rainfall rate (VRR) as a fmotion 
of T . . Solid curve calculated from Eq, [18) . 
Dotsare 30 minute average maximum VRR adapted 
from NA . 


Other factors affecting the precipitation 
efficiency include vertical shear and the environ- 
mental humidity. Empirically, efficiency decreases 
with increasing shear (Marwitz, 1972; Foote and 
Fankhauser, 1973). In our case study, the shear 
from cloud base to tropopause is 3.3xl0" 3 , which 
converts to an efficiency of approximately 35%, 
using the diagrams in the referenced papers. This 
is reasonably close to the average value over the 
range of our case study thunderstorms, but sig- 
nificantly below that for our intense storms. 

A dry environment decreases the efficiency 
in two main ways, as discussed by Brown and Klemp 
(1980) in their analysis of threerdimensional cloud 
model water budgets. Evaporation of the cloud edge 
and entrainment of the dry air reduces updraft size 
and intensity, and the cloud, therefore, processes 
less water. Also, the evaporation of rain before 
reaching the ground is greater with a drier environ- 
ment. Although the two shear conditions examined in 
Brown and Klemp's study were not indicative of the 
range of possible shears, their results indicated 
that the moisture effect on efficiency was more 
important than the effect of shear. 


3. INITIAL COMPARISONS WITH ONE-DIMENSIONAL 

MODEL RESULTS 

The comparison of one-dimensional model 
calculations with the simple formulations of 
section 3 and with the observations in NA 
has two general purposes. First, it is the initial 
step in trying to verify the simple relations in 
section 2 with more complicated convection models. 
Ultimately, the moisture budgets of simulated three- 
dimensional model storms must be examined to help 
determine the role of various factors (e.g., shear, 
environmental humidity). This work has been started 
by Brown and Klemp (1980). 

The second purpose has to do with the 
transportability of the relations in section 2. 

Eq. (16) seems to produce an upper bound for the 
data (see Fig. 4) and its linear form approximates 
that of the small data set. Therefore, the rain- 
fall is approximately proportional to the vertical 
mass flux, as indicated by the area average w . 
Therefore, if we can make an estimate of a ax 
storm's vertical mass flux, we can in turn make 
an estimate of the upper bound in rainfall (assum- 
ing knowledge of the moisture source). However, 
from satellite observations, our most common 
measure of relative updraft intensity is T . . 

With the present data set, we determined tWJ n T . / 
mass flux (w ) relation empirically (Fig. 5). 
However, thiSrSame relation will not hold under 
other climates or even under other synoptic condi- 
tions. A T . of 210K in mid- latitudes with a 
tropopause temperature of 21 5K implies a greater 
vertical velocity than the same T . in the tropics 
with a colder (higher) tropopause. n The one- 
dimensional model may provide the basis for deter- 
mining the T • /mass flux relation applicable under 
particular conditions, and therefore, help make 
Eq. (16) transportable to different environmental 
conditions and usable with the satellite-observed 

T min* 

The Simpson-Wiggert (1969) model was run 
on the Monett, MO sounding for 2315 on April 25, 
1975, for comparison with the data set from NA. 

Five runs for different updraft radii were 
made to produce clouds growing to various heights. 

A summary of results is shown in Table 1. The 
peak vertical velocity ranges from 14 to 50 ms". 
This is a reasonable range of vertical velocities 
for moderate to very intense thunderstorms. For 
comparison with the satellite cloud top ascent 
rates, the model w is ? converted to an area 
average over 100 x km , by assuming a uniform 
peak velocity over the updraft area and zero out- 
side of that. This produces the w in Table 1, 
which has values closer to those or^he satellite 
estimates which have a range from 1.2 to 8.0 ms"'. 

The model rainfall rate is calculated 
by dividing the model precipitation production by 
the time required for parcel ascent from cloud 
base to cloud top, in this case study estimated 
by cloud depth divided by w /2. The model VRR 
is simply the rainfall ratewll tipi ied by updraft 
area. Figure 7 displays model VRR versus w 
(model) where the line is drawn through the^del 
calculated points (x's) and the dots are the data 
from NA. The model results are very similar to 
that shown in Fig. 4 and again indicate that the 
VRR is proportional to vertical mass flux under 
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TABLE 1. Summary of Results for 
One-dimensional Cloud Model Runs 




Updraft Radius (km) 

1 


0.5 

1.0 

2.0 

3.0 

4.0 

Cloud 

Depth (km) 

8.6 

11.1 

12.6 

13.4 

13.8 

T »tn "<> 

217 

200 

188 

185 

182 

"max < mr1 > 

14 

22 

36 

45 

50 

"max (rnrl > 

0.1 

0.7 

4.5 

12.7 

25.1 

Rainfall , 
Rate (mm h" ) 

153 

340 

611 

774 

851 

35-i 

m s 

33 

297 

2133. 

6079. 

11882. 



Figure 7 . Volume rainfall rate (VFR) computed 
from the one- dimensional cloud mod^l output as 
a function of area-average (100 km ) model 
maximum vertical velocity (solid line connecting 
calculation points, x’s). Dots are radar VFF's 
and satellite cloud top ascent rates as in Fig. 4. 


the same environmental conditions. It will be 
interesting to determine if a similar simple 
relation exists with three-dimensional model 
output . 

The critical vertical veloci ty/T m1 * n 
relation is shown in Fig. 8. This curve is 
similar to that in Fig. 5 (satellite observations), 
except for the obvious difference that the model 
T . * s are significantly lower than the satellite 
ooservations. This is because the model T - n 
represents a very small horizontal dimensiAn at 
the center of cloud top, while the satellite T mi - n 
represents an average area over a large 
(100 km 2 ) area, therby being necessarily much 
warmer. In addition, there may be differences 
because the model T . is an updraft core tempera* 
ture while an infrared radiometer will measure 



Figure 8 . Model area-average w as a function of 
model T . and model T . (whicff 1 takes into 
account m k n temperature sfvtft to make the model 
output more comparable to the satellite observa- 
tions) . Data points are satellite observations 
from Fig . 5. 


cloud surface temperature. The relation of 
satellite T . to smaller scale T^ n has been 
examined by Adler et al. (1981) using aircraft 
overflight data with a radiometer similar to that 
on SMS/GOES. The results indicate that there is 
roughly a 14K difference between the satellite 
T . and the small scale (100 m x 100 m) T . due 
tS'the area average effect. Thus, a second hori- 
zontal scale has been added to Fig. 10 with 

T m1n (nradel) = T min (model) + 14 (1<0 

to allow a better comparison of the model results 
with the satellite observations (plotted as dots 
in the figure). The magnitude of the offset is 
based on only a few cloud observations on one case 
study day, and should not be considered general. 
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The NA satellite data set had T . 
ranging from 231 K to 201 K. This gives approximate 
agreement with the model results ranging from 
updraft radius 0.5 km to a radius between 2 and 3 km. 
The model results for w (even with the shift of T . 
axis) are low compared to the satellite cloud top min 
ascent rates (Fig. 8). The difference is very sig- 
nificant at the warm or less intense end of the 
range. The difference is unexplained at this time, 
but may be the result of the assumptions being made 
in order to make the variables comparable. An 

additional shift of 3K would produce significantly 
better agreement and is not out of the question 
according to Adler et al . (1981). However, the shape 
of the curve represents the trend of the data well. 

Figure 9 shows the curve representing 
the model-computed ViRR and the NA data set plotted 
on the T . scale. The model calculations handled 
in this way grossly underestimate the rainfall in 
the intense storms, but the curve has the sharp 
change in slope apparent in the data. 

Finally, Fig. 10 shows the VRR/T . 
relation when Eq. (16) is used and w is min 
substituted for by the w /T . mod?}f x relation of 
Fig. 9 instead of the sateTl ixeempirical relation. 
Again, the rainfall is underestimated, but not as 
badly. However, since the VRR/w relations of 
both Figs. 4 and 8 are more reasBftable, the failure 
must be in the w /T . connection. Obviously, 
more work must girfnto 1 Understanding the model 
results (both one-dimensional and three-dimensional) 
and their relation to the satellite observations. 



230 220 210 200 


T min (model) IK) 

Figure 9. Model volume rainfall rate (VRR) as a 

function of model T . urith the data points 

plotted on the T . scale . 

rmn 



235 230 220 210 200 

T MIN <K> 

Figure 10, Volume rainfall rate (VRR) calculated 
from Eq. (16) and using the model 
relation to convert to T . . Data points are 
same as in Figs. 6 and 9. 


A primary test of the use of the model 
results will be the testing of the relations under 
different climatological and synoptic conditions. 

One such test will be with radar results of maxi- 
mum echo height (convertible to T . ) versus rain- 
fall by Konrad (1977), Wilk (1981 T and Gagin and 
Lopez (personal communication) for various clima- 
tological regimes. Analysis by one of the authors 
(Atlas) already Indicates significant climatological 
differences. 


4. SUMMARY AND CONCLUSIONS 

The production of precipitation in 
thunderstorms is a complex confcination of cloud 
dynamics and microphysics and the interaction of 
the storm with its environment. However, in 
section 2 the simple steady state condensation 
rate equation is used to derive some formulations 
relating to volume rainfall rate (VRR) to satellite 
observables indicative only of relative updraft 
intensity. These relations reproduce the case 
study results of Negri and Adler (1981) to a first 
approximation. They emphasize that the VRR is 
roughly proportional to the area-average vertical 
velocity or mass flux (assuming the moisture source 
is constant). A satellite-based estimate of 
vertical mass flux in a thunderstorm, coupled with 
knowledge of the low-level moisture source (from 
surface conventional data) may be combined into an 
estimate of the total condensation rate--an upper 
bound for the VRR. The parameter T m i n can be used 
as a proxy for the vertical mass flux through an 
empirical relation. Thus, we have a physically- 
based relation between T m i n and storm rainfall. 
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Analysis of one-dimensional model results 
also indicates that model -based VRR is linearly 
related to the vertical mass flux. An important 
model result is that the model -generated relation 
between w and T m i n is similar to the empirical 
(satellite-based) relation and gives support to 
the idea of using the model results to establish 
Tmin/fna ss flux relations for various climatic or 
synoptic conditions, thus making the approach 
transportable. The idea of using the one- 
dimensional model to aid in interpreting the 
satellite observations appears plausible, but 
needs a deeper examination with reference to 
more sophisticated (three-dimensional) models 
and additional data sets. 

The authors have been pleased with 
the generally positive results presented in this 
paper. However, they are acutely aware that the 
simple relations tested here against one data set 
must be evaluated against additional sets for 
confirmation. Analysis of past radar data sets 
points out differences in rainfall /maxi mum height 
relations depending on climatological or, possibly, 
storm type variables. These will have to be taken 
into account for a successful result. 

Because of the complexity of the rain- 
fall process, it is unlikely that the approach 
presented in this paper will result in a deter- 
ministic satel 1 i te/rainfall relation. However, 
the approach may produce a basis for statistical 
relations for each climatological regime or 
storm type. 
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ABSTRACT 

Short-interval geosynchronous infrared satellite data are used to examine 1 1 cases of tornadic thunder- 
storms with respect to cloud-top temperature (height) variations relative to tornado touchdown times, 
and in three cases relative to the initial observation of mesocyclones by Doppler radar. The scale of the 
updrafts observable with the satellite infrared data is - 10 km. The cases are limited to those with relatively 
intense tornadoes. In 8 of the 1 1 cases there is a period of rapid cloud-top ascent 30-45 min prior to tornado 
touchdown. This upward growth appears to be associated with the formation of the mesocy clone. This 
ascent is followed by a period of no growth or even a drop in cloud-top height preceding, or at the time of, 
tornado touchdown. In the three remaining cases cloud-top ascent is evident in the satellite data at tornado 
touchdown. 


1. Introduction 

The study of thunderstorm top height changes 
near the time of tornadoes is important in order to 
help understand the relation of cloud dynamics to 
mesocyclogenesis and tornado touchdown. Also, 
since the storm top is the only portion of the cloud 
observable with satellite data, it is important to de- 
termine what tornado-related characteristics can be 
inferred from current satellite observations. 

Radar observations (e.g., Lemon et al . , 1978) 
indicate that tornado touchdown is often accom- 
panied by a decrease in echo maximum height and 
a decrease in the height of the Bounded Weak Echo 
Region (BWER). Both decreases are indicative of a 
weakening of the updraft. Fujita ( 1973) 2 proposed 
that the collapse of overshooting thunderstorm tops 
is associated with tornado generation and touch- 
down. Purdom ( 197 1) 3 showed examples of a pause 
in the rate of anvil expansion (measured from satel- 
lite image data) at the times of tornado touchdown. 
More recently, five tornado bearing clouds on 6 May 
1975 were discussed by Adler and Fenn (1979a). An 
examination of the five clouds elements, having 
eight tornadoes clearly associated with them, in- 
dicated that in seven of the eight cases the first 
report of the tornado took place during, or just after. 


1 Present affiliation: Office of Hydrology, National Weather 
Service, Silver Spring, MD 20010. 

2 Fujita, T.T., 1973: Proposed mechanism of tornado formation 
from rotating thunderstorms. Preprints 8th Conf. Severe Local 
Storms , Denver, Amer. Meteor. Soc., 191-196. 

3 Purdom, J. F. W., 1971: Satellite imagery and severe weather 
warnings. Preprints 7th Conf. Severe Local Storms , Kansas City, 
Amer. Meteor. Soc., 120-137. 


a rapid expansion of cloud-top areas defined by 
equivalent blackbody temperature ( T BB ) isotherms. 
This rapid expansion of areas within isotherms im- 
plies ascent on a scale of ~30 km, even larger than 
the scale of an individual infrared (IR) data point 
(10 km). 

In this paper, eleven tornadic storms, including 
five cases from Adler and Fenn (1979a), are evalu- 
ated with respect to cloud top temperature changes 
relative to tornado touchdown, and in three cases 
relative to initial reports of Doppler radar observed 
mesocyclones. 

2. Analysis of 11 tornadic thunderstorms 

Digital image data from SMS/GOES geosyn- 
chronous satellites were analyzed on four case study 
days (see Table 1) with the aid of an interactive 
computer/image analysis system called AOIPS 
(Atmospheric and Oceanographic Information Proc- 
essing System). Data were from times when the 
satellite was operated in the short-interval mode, 
images produced every 3-7.5 min over a swath of 
limited north- south extent. This high time resolu- 
tion is necessary to observe rapidly changing con- 
vective storms. The digital infrared (IR) data which 
are the basis of this study have a spatial resolution 
of —10 km (depending on location relative to the 
subsatellite point) and an effective temperature 
resolution of 1-2 K. 

The set of tornadic thunderstorms was limited, 
for this study, to storms with relatively strong tor- 
nadoes. All but one of the storms had tornadoes of 
intensity F2, or greater [on an intensity scale of F0 
through F5 as described by Fujita (1973)]. Kelley 
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Table 1 . Tornadic thunderstorm case study days. 


Date 

Area 

Spacecraft 

6 May 1975 

Nebraska 

SMS-2 

24 April 1975 

Missouri 

SMS-2 

20 March 1976 

Illinois 

GOES-1 

20 May 1977 

Oklahoma 

GOES-1 


et ai (1978) have shown from climatology that while 
only 38% of tornadoes are F2 or above, they ac- 
count for 98% of tornado fatalities. In a few cases 
where more than one tornado was produced from a 
storm, only the data relative to the time of the most 
intense tornado was examined. One storm with only 
an FI tornado was included in the sample because 
of its long tornado path. This procedure was de- 
signed to concentrate the analysis on so-called 
“supercell” storms (Browning and Foote, 1976), 
which are the main producers of relatively strong 


tornadoes (Lemon and Doswell, 1979). It also elimi- 
nates the relatively weak tornadoes, for which the 
time of occurrence is often suspect. 

The results of the 11 storms are shown in Fig. 1. 
If there is more than one tornado associated with 
a cloud, only the most intense tornado is considered. 
Such is the case in Figs, la, lc, Id and If. 

Fig. 1 depicts the variation with time of the satel- 
lite-observed minimum T BB for each storm relative to 
the time of the associated tornado. Decreases in T HB 
imply cloud-top ascent and increases in T BB imply 
descent. Because of the relatively coarse spatial and 
temperature resolution in the IR data, small or tem- 
porary changes in T BH or in the slope of the T BB curve 
should be considered suspect and may be artifacts of 
the data. For example, in Fig. li there is a brief stair- 
step at -40 min, which is probably a data artifact 
due to the coarse temperature resolution. These 
stairsteps are common in the data set and usually do 
not represent pauses in the ascent (or descent), as 



-60 40 -20 0 «20 -60 -40 -20 0 *20 -60 40 -20 0 *20 


TIME (MIN) TIME (MIN) TIME (MIN) 

Fig. 1 . Minimum cloud-top temperature (T HH ) as a function of time relative to the time of tornado touchdown in 1 1 cases of relatively 
strong tornadoes. The circled dots indicate the cases of the lowest temperature not being associated with a closed T Hi{ isotherm. The 
time of the tornado touchdown is noted by the vertical dashed line at the time equal to zero. In three cases (ij, k ) the initial times of 
mesocyclones observed by Doppler radar of the National Severe Storms Laboratory in Norman, Oklahoma are shown by the vertical 
dash-dot lines. 
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can be seen by the continued expansion of the num- 
ber of data points within the minimum T BB iso- 
therm during a period of constant minimum T BB . The 
relation of T BB isotherm expansion near cloud-top 
center to vertical velocity and divergence has been 
discussed by Adler and Fenn (1979a). 

Another problem is related to the coarse spatial 
resolution ( — 10 km) and the position of the satellite 
IFOV relative to the cloud top. A thunderstorm top 
will produce a slightly lower observed T IiB if the top 
is centered in an IFOV instead of being split between 
two adjacent views. Successive satellite views of a 
constant height storm may occasionally show an 
observed change in T BBi implying a false height 
change. An example of this may be the oscillation 
in panel / at - 10 min. 

In general, because of these data problems, small 
or short-lived changes in T HIi or its slope should be 
ignored. Changes of 3 K (4 K for the 20 May 1977 
storms) can be considered significant, along with 
smaller variations if they are supported by more than 
one data point at each end of the change. 

Eight of the 11 storms (Figs, la, lb, lc, Id, le, li, 
lj and lk) had an associated rapid decrease of cloud 
top T bb (as observed from satellite) -30-45 min be- 
fore tornado touchdown. This implies a rapid up- 
ward movement of the cloud top. These storms also 
had their period of most rapid area expansion within 
T bb isotherms prior to the touchdown. On most of 
these eight storms the rapid decrease in temperature 
was followed by a reduction in the slope of the 
temperature-time profile, and in the 20 May cases 
(Figs, li, lj and lk) a warming, indicating a cessa- 
tion in upward growth or a slight drop in maximum 
cloud height. The Fort Cobb case (Fig. lg) exhibits 
a substantial 6 K warming or cloud-top collapse. 

The sequence of events at the cloud top can be 
pictured as being associated with storm evolution 
in the following way. A period of rapid upward 
growth at cloud top is associated with or possibly 
precedes the formation of the mesocyclone at mid- 
levels in the storm. Tornado touchdown is usually 
associated with a decreasing or near constant maxi- 
mum cloud top. This sequence of events at cloud top 
is consistent with radar studies (e.g., Lemon et al . , 
1978) and with hypotheses proposed by Fujita 
(1973) 2 which were based on aircraft observations. 

Three storms (Figs. If, lgand lh) showed decreas- 
ing cloud top T bb (increasing height) at the time of 
the tornado touchdown. The tornadoes from all 
three of these storms were intense with Fujita/ 
Pearson scale ratings (Fujita, 1973) of F4, F2 and F4 
for the storms in Figs. If, lg and lh, respectively. 

The Neosho, Missouri storm (Fig. 1) is the most 
obvious exception to the other eight storms. While 
the tornado was on the ground, the cloud-top mini- 
mum temperature dropped nearly 10 K [over an 
area of —100 km 2 , one instantaneous field-of-view 


(IFOV)], indicating strong ascent. It is possible that 
the strong growth was associated with a new cell, 
adjacent to the tornadic cell, but this does not seem 
to be the case from the satellite analysis and exami- 
nation of available radar data. The cloud top ascended 
during the first 25 min of the tornado until the time 
the short-interval data ended. 

The Hoyleton, Illinois tornado (Fig. lg) of 20 
March 1976 also was on the ground while the satel- 
lite-viewed cloud top was ascending. However, the 
Hoyleton storm did show rapidly increasing height 
from 40 to 10 min before tornado touchdown, until 
a brief height plateau was reached at 207 K, just 
before tornado touchdown. During that constant 
temperature period, the area within the 207 K T BB 
isotherm decreased from 200 to 100 km 2 , which is 
equivalent to a change from two IFOV’s to one 
IFOV. This area decrease at constant temperature 
implies cloud top collapse, or at least constant 
height. This sequence is similar to an occurrance 
noted with the Omaha tornado on 6 May 1975 (Fig. 
8, Adler and Fenn, 1979a). However, in the Hoyle- 
ton case, the storm resumed its upward climb, and 
continued to ascend during the lifetime of the 
tornado. 

Fig. lh shows the Sadorous storm with rapid 
ascent during tornado touchdown, followed by a 
period of constant height and then descent. The 
largest area within the 206 K isotherm occurred at 
the first observation, and then decreased, indicating 
that maximum storm height happened just after tor- 
nado touchdown. Fujita et al. (1976) have analyzed 
this storm and indicate that the tornado reached F4 
intensity —20 min after touchdown. Perhaps, in this 
case, the cloud-height decrease was associated with 
the intensification of the tornado vortex and not 
touchdown itself. 

Therefore, a majority of the storms analyzed had 
rapid ascent of cloud tops 30-45 min before tor- 
nado touchdown, persumably during formation of 
the mesocyclone (documented in three cases). A 
minority of the storms, however, had ascending 
cloud tops (at the scale of the satellite observation) 
at the time of tornado touchdown. 

The indicated characteristic cloud-top evolution 
of the majority of tornadic storms means that we 
have a first link between satellite observations and 
tornadic activity. It does not mean, however, that 
we now have a way to discriminate between tornadic 
and non-tornadic thunderstorms with satellite data. 
Intense, non-tornadic storms examined also have a 
period of rapid growth above the tropopause, fol- 
lowed by relatively steady heights. 

The magnitude of the vertical velocities are rea- 
sonable when the horizontal scale (10 km) on which 
they are applicable is considered. A typical rate of 
temperature ( T BB ) decrease in Fig. 1 is 0.4 K min" 1 . 
Using an appropriate lapse rate, the temperature 
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change rate can be converted to an equivalent verti- 
cal velocity or cloud-top ascent rate. Hasler and 
Adler (1980) 4 have compared satellite observed 
height (from satellite stereo analysis) and IR tem- 
peratures for thunderstorm tops. For overshooting 
tops (above the tropopause) they find a temperature- 
height relation of 2.5 K km" 1 , which was approxi- 
mately midway between ambient and adiabatic for 
the particular case. However, the horizontal scale of 
the two observations are different, with 10 km for 
the IR temperature and 1-3 km for the stereo- 
determined height, depending on the size and 
definition of the feature in the visible data. If the 
two measurements were made at the same scale, 
the calculated lapse rate would shift toward adia- 
batic. Therefore, using the 2.5 K km 1 lapse rate 
in Eq. (1) 


produces a vertical velocity w appropriate on a hori- 
zontal scale of 1-3 km, roughly equivalent to radar 
resolution at thunderstorm top. Using a 2.5 K km 1 
value for the lapse rate gives an equivalent vertical 
velocity of 2.7 m s~ l . This is nearly equal in magni- 
tude to the echo top rise rate above the tropopause 
indicated for the Union City Storm (Burgess and 
Lemon, 1976 5 , Fig. 5.6; Lemon et aL , 1978, Fig. 13). 
However, larger cloud top ascent rates (up to 8 m s ') 
have been calculated from satellite IR data for in- 
tense convection penetrating through the upper 
troposphere (Adler and Fenn, 1979b). 

In Fig. 1 only panels li, Ij and lk show a T BB in- 
crease (indicating descent or cloud top collapse) 
prior to, or at the time of, tornado touchdown. The 
6 May cases (Figs, la- lc) show a constant tempera- 
ture with time at the time of the tornado. This dif- 
ference (between the two days) may be related to the 
difference in larger scale cloud structure between 
the two cases and the satellite sensor's response 
characteristics. The 6 May storms appeared in the 
satellite images as a narrow line lying along a gener- 
ally north-south oriented cold front. The east- 
west extent of the storm's anvils increased with 
time. As the satellite IR sensor scans from left to 
right across the scene it moves abruptly from a 
warm target (ground) to a cold target (thunder- 
storm top). With a narrow cirrus shield, the sensor 
may not have enough time to accurately respond to 
the minimum T BB , especially since it is often located 
near the left, or upwind side of the anvil. That is, 


4 Hasler, A.F.,andR. F. Adler, 1980: Cloud-top structure of a 
tornadic thunderstorm from 3 min interval stereo satellite images 
compared with radar and other observations. Preprints 19th Conf. 
Radar Meteorology , Miami, Amer. Meteor. Soc., 405-412. 

' Burgess, D. W., and L. R. Lemon, 1976: Union City Storm 
history, the Union City, Oklahoma tornado of 24 May 1973, 
R. A. Brown, Ed. NOAA Tech. Memo ERL NSSL-80, 35-51. 


for narrow (in the left-to-right direction) cirrus 
anvils, the minimum T BB is probably overestimated 
(i.e., estimated to be warmer than it is). However, 
as the originally narrow anvil expands with time, 
the amplitude of the temperature overestimation 
should decrease as the distance from the left anvil 
edge to the storm center increases. If, while the anvil 
is expanding, the actual minimum T BB is warming 
slightly, this could be obscured by the varying over- 
estimation effect, with the result being a near-con- 
stant observed T BB with time. For example, cloud 4 
in Fig. la has a 35 min period of constant T BB (212 K) 
while the anvil is expanding. The distance from the 
left anvil edge (denoted by the 226 K isotherm) to 
the center of the area of coldest T BB increases from 
two data points to nine data points over the 35 min 
time period. Work on a technique to remove this 
sensor response effect is planned. 

On 20 May 1977, the storms (Figs, li, lj and lk) 
were imbedded in a very large cirrus shield with a 
large distance between the left edge and the storms’ 
centers. Thus, the sensor response problem in terms 
of estimating the minimum T BB should not exist. 
In these cases, a warming was observed. 

3. 20 May 1977 storms 

The three intense tornadic storms observed on 
20 May 1977 are the best analyzed storms to this 
date with respect to Doppler radar coverage. A com- 
plete analysis of the cases with satellite, Doppler 
radar reflectivities and velocities, and other data is 
underway by the authors and others. Preliminary 
results based primarily on satellite data will be dis- 
cussed in this section. 

a. Alt us storm 

The Altus, Oklahoma storm was actually a com- 
bination of two related thunderstorms which could 
be identified in both the satellite and radar data. The 
National Severe Storms Laboratory (NSSL) Doppler 
radar log indicated only one mesocyclone, but an 
examination of the positions of the mesocyclone, 
reflectivity maxima and satellite features showed 
that two mesocyclones, associated with two storms, 
were involved. The temperature-height versus time 
diagram for both cells is shown in Fig. 2. Cloud 37 
(the numbering system is from a larger study) was 
first defined as penetrating the cirrus shield just after 
1840 GMT (1140 CST) at 207 K. The background 
cirrus T BB was 209 K, approximately equal to the 
tropopause temperature. With dense cirrus shields 
already existing, such as in this case, the new 
thunderstorms cannot be detected in the satellite 
IR data until they significantly penetrate the cloud 
deck at very cold temperatures. 

Cloud 37 peaked at 198 K and remained well de- 
fined at 200 K as it traveled north-northeastward. 
A mesocyclone was first noted at 1932 GMT in the 
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Doppler radar observation log, colocated with the 
satellite feature. The radar-observed circulation, 
marked as M37 in Fig. 2, perhaps existed earlier, 
but was not detected because of the relatively long 
distance (220 km) from the radar to the storm. This 
mesocyclone remained relatively weak, and was dis- 
sipating at the 1952 GMT observation. 

At about 1940 (all times GMT unless indicated 
otherw ise), evidence of a second storm (Cloud 44) 
was observed in the satellite data approximately 
30 km to the southeast of cloud 37 and colocated 
with a second radar reflectivity maximum. Cloud 44 
ascended rapidly to 200 K (~4 km above the tropo- 
pause) by 1950. An associated mesocyclone was first 
observed at 1952 by the Doppler radar. Thus, the 
initiation of the cyclonic rotation appears associ- 
ated with rapid cloud top ascent and, therefore, a 
strong updraft. After a relatively stable 15 min 
period the cloud top temperature oscillated and then 
warmed. The oscillation around 2010 may not reflect 
an actual cloud height fluctuation, but an effect of 
sensor resolution and the location of the satellite 
data point relative the highest part of the cloud top. 
However, between 2005 and 2020 there was a slight 
storm top decrease, just preceding the reported tor- 
nado touchdown. The storm top ascended again just 
after tornado dissipation, a feature also seen in the 
following two cases. 

b . Fort Cobb storm 

The Fort Cobb, Oklahoma storm developed to the 
southeast of the Altus storm and also had two as- 
sociated cold features in the satellite IR data. One of 
the features, identified as cloud 49, was first defined 
at 2045 and reached a temperature of 202 K, then 
weakened and became difficult to define in the IR 
field until 2210. At that time (see Fig. 3), the time 
sequence of IR images showed a sharp cloud top 
ascent from 204 to 198 K in 18 min. A mesocyclone 
(mesocyclone A) was initially observed at the time 
of maximum height. The next observation of the 
storm was 12 min later and the storm top has warmed 
4 K, or ~2 km, based on a lapse rate halfway be- 
tween ambient and adiabatic. After 2240 the cold 



Fig. 2. Minimum cloud top temperature as a function of time 
for cloud elements associated with the Altus, Oklahoma tor- 
nado of 20 May 1977. 



Fig. 3. Minimum cloud top temperature as a function of time 
for cloud elements associated with Fort Cobb, Oklahoma storm 
of 20 May 1977. 


feature in the IR data became more difficult to follow 
and may have been just a piece of outflow debris. 
The times of the associated tornadoes are shown in 
Fig. 3 based on analysis presented by Johnson et al. 
(1980). 6 The touchdown of the Fort Cobb tornado 
occurred at the end of a period of sharp decrease in 
the satellite-observed height with a temperature 
change from 198 to 204 K. 

As the first IR feature (cloud 49) became indis- 
tinct and disassociated from the main radar echo 
(and mesocyclone A) at approximately 2230, another 
T bb minimum (cloud 49a) was defined 30- 40 km to 
the southwest of the mesocyclone position. This 
new feature in the IR field had an associated radar 
reflectivity maximum. Between 2230 and 2300, the 
satellite data indicated cloud-top ascent as this new 
updraft center (as defined in the T BB field) raced 
toward the northeast at 30 m s’ 1 , finally catching up 
to the location of the original reflectivity maximum 
at approximately 2250. The new IR feature, labeled 
49a, became dominant and was located near both 
mesocyclones (which are approximately 8 km apart), 
although it was probably related to the updraft as- 
sociated with mesocyclone B. A detailed study of the 
complex evolution of radar and satellite-observed 
features during this period is underway. 

The top of cloud 49a continued to ascent until 
2310. The apparent flattening of the curve between 
2255 and 2305 is an artifact of the poor vertical 
(temperature) resolution (2 K) in the satellite data. 
Although the minimum temperature remained con- 
stant at 202 K the number of IR data points at that 
temperature increased from two to ten, indicating 
continued ascent. The second tornado, associated 
with mesocyclone B and an updraft along the gust 
front (Ray et al . , 1981), touched down during this 
period of apparent cloud-top ascent. Again, as in the 
Altus case, the cloud top rose after tornado dis- 
sipation. 


6 Johnson, B. C., K. W, Johnson, P. S. Ray, J. S. Bradberry, 
J. J. Stephens and W. C. Bumgarner, 1980: The morphology of 
some tomadic storms. Preprints / 9th Conf. Radar Meteorology , 
Miami, Amer. Meteor. Soc., 311-316. 
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c. Del City storm 

The Del City storm also involved two thunder- 
storms, identified in the satellite data as clouds 53 
and 56 (see Fig. 4). Cloud 53, which had a period of 
rapid growth at 2240 to a minimum temperature of 
200 K is referred to as the “hailstorm’’ in the study 
by Ray et al. (1981). As cloud 53 collapsed another 
storm (cloud 56) appeared to the south, grew, moved 
northward, and caught up and merged with the hail- 
storm. Cloud 56, the main feature, reached a tem- 
porary height maximum at 202 K at 2305-2325, then 
descended slightly to 204 K. A new surge in the up- 
draft was noted between 2335 and 2355 as the storm 
top reached 200 K. This period of cloud top ascent 
was coincident with the formation of the meso- 
cyclone. 

The cloud-top temperature fell back to 202 K, 
then remained nearly constant for one hour as the 
storm moved in a northeast direction. Although this 
period was marked by constant temperature, the 
position of the cold feature was twice redefined 
slightly southeastward in a non-continuous manner. 
This redefinition is probably related to the effect of 
the IR cold area beginning to represent debris from 
the updraft and moving off from the storm core and 
radar echo until a new surge in the vertical velocity 
reestablishes the IR cold point over the updraft 
itself. 

No cloud-top collapse was noted in association 
with tornado touchdown in this case. Although there 
is a 30 min absence of data centered at 0015 GMT, 
there is evidence supporting cloud top rise after tor- 
nado dissipation. This was also evident in the Altus 
and Fort Cobb cases, and may represent reintensifi- 
cation of the storm through a new updraft on the 
storm flank, or revitalization of the original updraft. 

4. Mechanisms for rotation generation and intensifi- 
cation 

The mechanisms by which certain thunderstorms 
produce the initial storm-scale rotation, or meso- 
cyclone, and the mechanisms by which this rotation is 
concentrated to produce a tornado are not completely 


understood. However, recent research using obser- 
vations from Doppler radars (e.g., Heymsfield, 1978) 
and results of three-dimensional numerical models 
(e.g., Schlesinger, 1980) have indicated that the two 
dominant terms in the vorticity equation related to 
vorticity changes following a fluid parcel are the 
concentration or stretching term, and the tilting 
term, shown in Eq. (2) as the two terms to the right 
side of the equal sign: 


di 

dt 


-£(V„-V) + 


d\v du 
dy dz 


dw dv \ 
dx dz j 


( 2 ) 


The vertical component of relative vorticity is given 
by £, while //, v and w are the components of the 
vector velocity V. This formulation neglects the 
solenoidal term, the frictional torque term, and 
terms involving the earth's vorticity. Although the 
stretching and tilting terms are the most important, 
Lemon and Doswell (1979) suggest that the solenoi- 
dal term may also be a significant, but secondary, 
contributor. 

The satellite observations presented in this paper 
cannot, by themselves, answer the vorticity mecha- 
nism question, but do complement the studies based 
on radar observations and other data. A majority 
(8/11) of the satellite cases indicate that during the 
development of the mesocyclone the thunderstorm 
top is rapidly ascending, implying intensification 
of the storm updraft. This relation confirms radar 
case study observations (e.g.. Lemon et al . , 1978) 
of echo top height increases during mesocyclonesis. 
Lemon and Doswell (1979), in developing a descrip- 
tive model of mesocyclone evolution, indicate that 
the Bounded Weak Echo Region (BWER), which 
implies updraft, is colocated with the mesocyclone 
in its early stages. If the vorticity center is indeed 
colocated with the updraft at this stage, the implica- 
tion is that the vorticity increase associated with 
mesocyclone generation occurs through the vertical 
stretching of the air column in the rapidly intensify- 
ing updraft. However, the tilting term must also be 
important, at least along the updraft edges, where the 
gradient of vertical velocity is large. Although nu- 
merical model results (e.g., Schlesinger, 1980) in- 



Fig. 4. Minimum cloud top temperature as a function of time for cloud elements associated 
with the Del City, Oklahoma storm of 20 May 1977. 
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dicate dominance by the tilting term in this early 
stage, this dominance is artificial because the models 
eliminate stretching at the earliest time by assuming 
no ambient vorticity (vertical component) in the 
initial flow field. 

Tornado touchdown was associated (in a majority 
of cases) with satellite-observed constant or de- 
creasing cloud top height. These observations are 
in agreement with radar echo-top information such 
as in the Union City tornado case (Lemon et al . , 

1978) . Lemon and Doswell (1979), synthesizing 
observations from a number of cases, conclude that 
the vortex intensification and tornadogenesis is 
associated with updraft weakening, downdraft for- 
mation upwind of the updraft and a shift of the vortex 
center from the updraft center to the updraft-down- 
draft boundary. This implies that the tilting term may 
be dominant at this stage, as weakly indicated by 
observations (Ray et al . , 1976; Lemon and Doswell, 

1979) , and strongly indicated in model results 
(Schlesinger, 1980). The satellite data merely em- 
phasize that vortex intensification is apparently 
related to updraft weakening in most cases. Al- 
though this observation also fits the “ vortex valve” 
concept of vortex-updraft interaction (Lemon, 
1976, 7 Lemon *7 a/. , 1978; Brandes, 1978), observa- 
tions of the intense vortex being situated on an 
updraft-downdraft boundary (e.g., Brandes, 1978) 
appear to eliminate this theory. 

The majority of storms depicted in Fig. 1, there- 
fore, agree with previous observations of echo top 
behavior in tomadic storms and the conceptual 
model of Lemon and Doswell (1979). The three ex- 
ceptions (Figs. If, lgand lh) may be just slight varia- 
tions of the model (this is plausible with Figs, lg 
and lh), or may be radical departures (likely for 
Fig. If). The Neosho term (Fig. If) displays rapid 
cloud-top ascent during and following touchdown 
of a very intense tornado (F4). This may imply that 
the stretching term may still be very important at this 
stage, in this case. 

5. Summary 

Infrared geosynchronous satellite data have been 
analyzed to determine thunderstorm top height vari- 
ations relative to the time of tornado touchdown. 
In eight of eleven cases of relatively strong tor- 
nadoes, there was a period of rapid height increase 
(temperature decrease) 30-45 min prior to tornado 
touchdown. A typical value for this temperature 
decrease above the tropopause is 0.4 K min -1 , which 
is equivalent to a cloud top ascent rate of ~3 m s' 1 . 
In three cases where Doppler radar observations 


7 Lemon, L. R., 1976: Tornadic storm evolution: Vortex value 
hypothesis. Appendix G, the Union City, Oklahoma Tornado of 
24 May 1973, R. A. Brown, Ed. NOAATech. Memo. ERL NSSL- 
80, 229-234. 


were available, this period of rapid cloud top ascent 
coincided with the development of the mesocyclone. 
In the same eight cases, the time of tornado touch- 
down was a period of nearly constant height, or 
slight height decrease (temperature increase). This 
sequence of satellite-observed events in the evolu- 
tion of these tornadic thunderstorms is similar to 
radar-observed echo top changes noted in other 
tornadic thunderstorms and is consistent with a con- 
ceptual model (Lemon and Doswell, 1979) in which 
the mesocyclone vorticity is generated primarily 
through the stretching of the vertical air column in a 
rapidly intensifying updraft. The intensification and 
descent of the vortex to form a tornado is related 
to a weakening of the updraft, formation of a down- 
draft, and a shift of the vortex to the updraft-down- 
draft boundary, so that the tilting term becomes 
dominant in the generation of vorticity. 
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OBJECTIVES 

Investigate the cloud top height structure of the 
Wichita Falls tornadic storm of April 10, 1979 
(and other tornadic storms on this day) from 0.5 
micron observations from two satellites (stereogra- 
phy). 

Identify individual convective elements in the 
10 km resolution, 3 min interval infrared data and 
monitor the temporal change in minimum cloud-top 
temperature (CTT) , with particular regard to time 
of tornado genesis and echo area expansion. 




Robert Mack 

and GE/MATSCO 

Beltsville, MD 20705 


Investigate the limitations imposed by infrared 
field of view in defining individual thunderstorms 
in GOES digital imagery. 

Compare stereographically determined cloud-top 
heights with simultaneous infrared observations and 
low-level reflectivity data to better understand 
the dynamics of storm development . 


Figure 7. Cloud- top temperature (CTT) analyses 
(K) of storms on April 10, 1979, 2252-2325 GMT. 
Areas of local temperature maxima and minima are. 
identified by crosses, while interior temperature 
maxima (warm regions) are shaded . Analysis of 
individual thunderstorms often requires 3 min 
interval data and full radiometric resolution (IK). 
The CTT minima associated with the Wichita Falls , 
Pilot Point and Harrold tornadoes are respectively % 
the points 212K, 212K, and 206K , aligned southwest 
to northeast in +he 2252 GMT panel. 


Figure 2. Same analysis as Figure 1 excert 
2334-0016 GMT 10-11 April 1979. A pronounced 
region of warmer temperature (21 7- 22 IK) has formed 
downwind of a region of colder (208-213K) "V" 
shaped region. The minimum temperature of this 
region is the storm top associated with the 
Seymour -Wichita Falls tornado. Possible explora- 
tions for this warm region are: 

a) an eddy produced by the obstruction of 
the environmental flow around the main updraft , 
entraining warmer cloud material from the anvil 
edge, 

b) warming produced by subsidence of air in 
the lee of the ascending tower , and 

c) ejection of cirrus particles by the over- 
shooting top into the ( warmer ) stratosphere . 
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Figure 3. Comparison of minimum CTT and radar 
reflectivity area for the Crowell -Lawton tomadic 
storm. Coincident with the appearance and 
expansion of the level 4 (45 dBz) echo area was 
a cooling of the minimum CTT from 21 IK to 206K in 
3 min . The first reported, tornado (at Crowell) 
was preceded by a slow cooling (4K/27min) , 
followed by a 2K warming. After 2230 GMT, both 
the level 4 echo area and the minimum CTT became 
indistinguishable from similar areas of neighbor- 
ing storms. Occasionally a temperature minima 
will dissappear into the broader temperature 
structure of the anvil. These points are denoted 
by an asterisk (*) . Tropopause temperature from 
OKC at 2000 GMT is 21 5K . 



Figure 4. Comparison of minimum CTT and radar 
reflectivity area for the Harr old- Grand field 
storm. A significant cooling of the cloud- top 
(8K/27 min) was followed by a 2K warming gust 
prior to tornado onset. This was followed by 
another period of cooling (5K/24 min) during which 
the level 4 (45 dBz) echo area increased by an 
order of magnitude and merged with another storm. 
Pulsations in the level 5 (55 dBz) echo area after 
2230 GMT appear to be manifested at cloud top y . 
with distinct episodes of cooling (1-3K) associa- 
ted with each maximum in echo area. Small tempera- 
ture changes at cloud tpV> uhen considered with 
respect to the large (100 km 2 ) view of the IP 
sensor 3 may be indicative of substantial changes 
in storm volumetric rain rate as observed by radar 
Minimum CTT for this storm was 203K 3 coincident 
with the largest areal extent of the level 5 
echo at 2341 GMT. 
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Figure 5. Comparison of minimum CTT and radar 
reflectivity area for the Seymour-Wi chita Falls 
tomadic storm. A cloud-top ascent (cooling) 
rate of 9K/46min followed by a warming of 2K/0 min 
preceded the onset of the Seymour tornado. A 
cooling of 7K/21 min followed, by a 3K/3 mi,n 
warming preceded the Wichita Falls tornado . The 
rapid fluctuation in minimum CTT between 2300 
and 0000 GMT appears to be related to the JR 
field of view viewing a small feature on the 
extreme upwind anvil edge. This would lead to an 
overestimate (warming) of the minimum CTT. (See 
discussion of Figures 7-10.) 
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Figure 6. Comparison of minimum CTT and radar 
reflectivity area for the Pilot Point tornadio 
storm, Thie storm underwent two periods of cooling 
and warming , one prior to, the other after, the 
identification of a hook echo in the low-level 
reflectivity . This was apparently a short-lived 
and weak tomadic storm . 


Figure 8. Comparison of infrared derived tempera- 
ture (K), stereo derived height (km) and low-level 
radar reflectivity for the cross-section D-D 1 in 
Figure 7 ( 2345 GMT) . Distance scale starts at the 
northern anvil edge . These mature, downwind storms 
reveal a narrow, intense band of precipitation 
centered beneath a symmetric region of CTT /height. 
Nearly colocated peaks are 55 dBz, 209K , and 
12.7 km . 



Figure 7. Stereo graphically derived cloud- top heights, (km, MSL) from GOES East and West images at 2345 
GMT on April 10, 1979 . The satellites were inadvertently synchronized to within 4s. Highest top on the 
image was that associated with the Wichita Falls storm (14.6 km). Absolute accuracy of the heights is 
approximately 1 km, relative accuracy about 0.5 km. Cross -sections of height, temperature and low-level 
radar reflectivity are portrayed in Figure 8 (Downwind, D-D ’), Figure 9 ( Upwind , V-U 1 ) and Figure 10 
(SW-NE section). 
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Figure 9 . Similar analysis as Figure 8 for the 
cross-section U-U* in Figure 7. Maximum reflec- 
tivity in this section was 35 dBz 3 located at the 
midpoint of a broad area of constant cloud- top 
height (13 km ) , but also located beneath the local 
warm area (220K) which formed in the lee of *he 
" V-Notch . " The data implies that above the . 
tropopause (21SK) no unique relationship exists 
between cloud-top height and temperature . 



Figure 10. Similar analysis as Figure 8 for the 
cross-section SW-NE in Figure 7. The Wichita 
Falls 3 torm (SO km) showed excellent spatial 
agreement between maximum reflectivity (55 dBz ) , 
maximum height (14.8 km) and minimum CTT (212K ) , 
however y significant rain extends into the warm 
region (219K). The coldest CTT at this time was 
20 7K for the storm at 250 km y associated with a 
local height maximum 3 but not with the absolute 
height maximum. It is postulated that the IE 
sensor is underestimating the magnitude of the 
convection (warm bias) for storms at the extreme 
upwind anvil edge. The data suggests that while 
points of maximum rainfall may be inferred from 
cloud-top height maxima (CTT minima) no unique 
relationship exists between all pixels at a given 
height (temperature) and rainfall beneath that 
pixel. 


CONCLUSIONS 

Analysis of individual thunderstorms in GOES IR 
imagery for severe weather applications requires 
3 min interval data at full radiometric resolution 
(one degree at cold end) . 

The farthest upwind storm (Wichita Falls storm) was 
denoted by a minimum in temperature located at the 
apex of a "V" shaped cold area. A local warm area 
developed immediately downwind of this tower, 
ultimately attaining a temperature 12° warmer than 
the ascending tower at 2346 GMT. 

Three explanations were proposed for the downwind 
warm area: an eddy formed in the wake of air 

diverging around the ascending tower; subsidence 
of air in the lee of the tower; and ejection of 
cirrus particles into the stratosphere. 

Immediately preceding each of five distinct 
tornado touchdowns was a rapid cooling, then 
warming, of the minimum cloud- top temperature. 

Small temperature changes in CTT above the tropo- 
pause may be indicative of fluctuation in storm 
volumetric rain rate as observed by radar. 


Stereographic analysis revealed the Wichita Falls 
storm top to be at 14.6 jjO.5 km at 2345 GMT, the 
highest point on the image. Minimum CTT at this 
time was 212K, 7° warmer than the coldest CTT; an 
apparent overestimate of temperature for a small 
feature near the upwind anvil edge. 

No stereo height minimum could be found in the 
downwind warm area, implying no unique relation- 
ship between ( temperature/height above the tropopause. 

Cross-sectional analysis suggests that while areas 
of maximum rainfall may be inferred from cloud-top 
height maxima, no relationship exists between all 
pixels at a given height and rainfall beneath that 
pixel. 
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ABSTRACT 

A relationship between maximum winds and satellite-measured equivalent blackbody temperatures 
near tropical cyclones is investigated with data from both the Atlantic and western North Pacific areas. 
This investigation revealed not only a significant correlation between satellite-derived equivalent 
blackbody temperatures and maximum winds but also a strong lag relationship between these 
temperatures and maximum winds. From this latter relationship a regression technique was developed 
to forecast 24 h changes of the maximum winds for weak (maximum winds *£ 65 kt) and strong (maximum 
winds > 65 kt) tropical cyclones by utilizing the equivalent blackbody temperatures around the 
storm alone, and together with changes in maximum winds during the preceding 24 h and the current 
maximum winds. Testing of these equations with independent data showed that the mean errors of fore- 
casts made by the equations are lower than the errors in forecasts made by persistence techniques. 


1. Introduction 

Hurricane caused damages in the United States 
average over $600 million per year (Gentry, 1966). 2 
Damages from tropical cyclones are even much 
greater in many other countries on the western 
borders of tropical cyclone basins when expressed 
as a percentage of the gross national product. 
Forecasts of where and when a hurricane will 
strike are obviously of great importance. 

Observations and forecasts of the maximum winds 
in hurricanes are also very important because dam- 
ages caused by hurricanes vary exponentially with 
the maximum wind speeds. While the force of the 
wind varies with the square of the speed, some 
of the historical surveys of total storm damage sug- 
gest that the damage varies with a higher power of 
the wind speed, i.e., 


1 Formerly associated with General Electric Space Division 
(MATSCO), Beltsville, MD 20771. 

2 Updated using records at the National Hurricane Center, 
Miami, Florida. 


D = KV\ (1) 

where D (total damage) is the total damage caused 
by the storm, K a constant, V the maximum wind 
speed, and n some number between 2 and 5 (Howard 
et al . , 1972). This relationship emphasizes the 
importance of knowing and predicting the intensity 
of a tropical cyclone. Furthermore, contrary to the 
improvements made in predicting tropical cyclone 
movement, relatively little progress has been made 
in the last 20 years in developing improved ob- 
jective techniques for forecasting maximum winds 
of these storms. Except in special situations, fore- 
casters still rely very heavily on persistence and 
climatological techniques when forecasting tropical 
cyclone intensity. As a result, operational skill in 
forecasting maximum wind speeds in hurricanes 
remains low. 

In spite of the need for knowledge of tropical 
cyclone intensity by the hurricane forecast services, 
aircraft reconnaissance of these storms on the aver- 
age is being reduced for economic reasons. Efforts 
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have been increased in recent years, therefore, to 
use satellite data to observe and predict the intensity 
of tropical cyclones. Results from these efforts have 
been very encouraging and they keep improving as 
satellite data of better quality become more readily 
available with each new satellite series. 

While satellite and other platforms usually pro- 
vide sufficient information to identify the current 
intensity of tropical cyclones, relatively few satellite 
data (visible, infrared or passive microwave) have 
been utilized in predicting changes in intensity, and 
then primarily in a qualitative sense. Satellite 
measured equivalent blackbody temperatures ( T BB ) 
of cloud tops in tropical cyclones are believed to 
contain useful information about storm intensity and 
expected changes of intensity. Latent heat released 
when the warm moist tropical air ascends in major 
cumulus towers of hurricanes is the primary fuel 
for the storms (Dunn and Miller, 1960), and its 
availability is indicated by the amount and vigor of 
the convection within the cyclone. This can be de- 
duced from satellite T BB measurements. 

The hurricane is a prolific producer of clouds. 
The convective towers build far into the troposphere 
and sometimes penetrate the lower stratosphere, 
thus producing very cold cloud tops. The high-level 
shearing and outward spiraling winds spread the cold 
cirrus over a large area beyond the region of most 
active convection. This air subsides as its spirals 
away from the storm center causing the cirrus to 
begin dissipating as the air warms adiabatically. 
These effects are easily observable in satellite 
imagery and can be quantified through measure- 
ment of cloud-top temperatures. Thus, areal 
distribution of T BB provides information on the 
extent and strength of the convection which serve 
as indices of the latent heat released and indicate 
the extent to which the clouds of the storm are 
organized into patterns. 

The latent heat is ultimately converted to the 
kinetic energy which causes the extreme winds of 
the tropical cyclones (Riehl, 1954). For this to take 
place, however, complex processes are involved in- 
cluding, among other things, conversion of the heat 
to potential and available energy. Finally, the 
kinetic energy has to be concentrated by the flow 
patterns usually into relatively narrow bands, for 
the storm to become truly destructive. All these 
processes take time and there should be a lag 
between changes in convective activity and changes 
of maximum winds in the storm. 

The results of theoretical-numerical model experi- 
ments simulating development and maintenance of 
tropical cyclones support the reasoning and suggest 
that maximum vertical motion, i.e., maximum con- 
vection precedes the highest winds by 1-3 days 
(Rosenthal, 1978; Kurihara and Tuleya, 1974). 
Riehl (1954) and Rosenthal (1978) have also empha- 


sized that the convection needs to be organized 
by some larger scale system into a suitable pattern 
(e.g., spiral bands and eye- wall) before rapid 
intensification of the tropical cyclone takes place. 

Dvorak and earlier investigators at the National 
Environmental Satellite Service have developed 
techniques to use satellite imagery to identify the 
present intensity of the tropical cyclone and to sug- 
gest future changes of the intensity (Dvorak, 1975; 
Hubert et al., 1969). While these techniques have 
shown skill and the latest Dvorak technique is in 
widespread use in the tropics worldwide, it still 
involves considerable subjectivity especially in the 
forecasting of storm intensity. Dvorak (1975) utilized 
the degree of pattern organization to identify the 
current storm intensity from satellite imagery. He 
found that the size of the central dense overcast 
of cirrus and the degree to which the spiral bands 
of convective clouds encircled the storm center to be 
important factors. Others who recently have studied 
the relationship between satellite measurements of 
clouds and tropical cyclone intensity, convection 
or rainfall include Arnold (1977) and Griffith et 
al. (1978). 

Based on the heuristic reasoning just presented, 
results from the theoretical experiments, results of 
using Dvorak's technique under operational condi- 
tions, and other research, the authors have de- 
veloped a hypothesis tested by the experiments re- 
ported in this paper. It says 1) the T BB of the tropical 
cyclone cloud tops provide a measure of the con- 
vection and an index of the latent heat released 
for eventual conversion into kinetic energy; 2) the 
T bb areal distribution serves as an index of the 
organization of the storm’s convective activity; and 
3) the lower (higher) the mean T BB of the cloud tops 
over a moderate sized area, the stronger (weaker) 
and more (less) persistent is the convection and 
the more likely that the maximum winds in the storm 
will increase (decrease) with time. 

2. The data and the analysis 

The infrared(11.5 /xm, window channel) T BB ior a 
number of tropical cyclones over open ocean areas 
were analyzed using data from the western Atlantic 
for 1969 and the western Pacific for 1970, 1973 and 
1974. The 1969 T BB were measured by the Medium 
Resolution Infrared Radiometer (MRIR) sensor on 
Nimbus 3 with a spatial resolution of 55 km at nadir. 
Those in 1970, 1973 and 1974 were measured by 
the Temperature Humidity Infrared Radiometer 
(THIR) sensor on Nimbus 4 and 5 with a spatial 
resolution of 8 km at the subpoint. Further details 
of the Nimbus 3, 4 and 5 satellites and their instru- 
mentations are provided in the Nimbus 3, 4 and 5 
User's Guides (Nimbus Project 1969, 1970, 1972). 
These data were analyzed using the scheme illustrated 
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Fig. 1. Grid used in analyzing the temperature data. The con- 
centric rings are spaced 111 km apart. The center of this grid 
coincides with the center of the tropical cyclone. 


in Fig. 1 in order to obtain a measure of the intensity, 
expanse and organization of the storm. The con- 
centric circles are 111 km apart and the rings they 
bound are numbered outward from 1 to 12. The 
mean T BB was computed for each ring with the 
center of the diagram coinciding with the center of 
the storm. In addition, the mean T BB was calculated 
for each octant of each ring (hereafter referred to 


as a sector) with the top of the diagram being 
oriented both toward the north and also along the 
direction of motion of the storm. 3 To get a further 
measure of how well the convective towers were 
distributed symmetrically and concentrated about 
the storm center, the standard deviations of the 
mean sector temperatures were computed for rings 
1-5 and for various combinations of rings. With 
these data it is feasible to study the expanse and also 
the organization of the storm as well as the intensity 
of the convection. 

The wind data for each Atlantic storm were ob- 
tained from the best track records maintained at 
the National Hurricane Center. The wind data for 
each Pacific storm were obtained from the best 
track and warning information compiled by the Joint 
Typhoon Warning Center, Guam (U.S. Fleet 
Weather Central, Joint Typhoon Warning Center, 
1970, 1973, 1974). A linear interpolation of wind was 
made when satellite observation time occurred 
between the best track data times. 

3. T bb and storm intensity 

The first tests made in this investigation were 
with 1969 Atlantic tropical cyclones. This investi- 
gation determined how well the cloud top T BB 
demonstrated an index of convection. The mean 
data for the rings composited for 16 hurricanes are 
compared with similar data for 19 storms of less than 
hurricane intensity (three of the tropical cyclones 
were still at depression stage) in Fig. 2. For rings 
1-4 (Fig. 1) the mean T BB were 7- 10°C lower in the 
hurricanes than in the weaker storms, but in rings 


3 Since all results reported in this paper used the entire ring of 
data, storm orientation is of no concern. 



Fig. 2. Comparison of mean temperatures of cloud tops around 16 hurricanes 
with mean temperatures from 19 storms of less than hurricane intensity, all in 
1969. The tropical cyclones all occurred in the western North Atlantic in 1969. 
The A T values are positive when the temperatures in the hurricanes are higher. 
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Fig. 3. Comparison of mean temperatures of cloud tops around 
1970 tropical cyclones of the western North Pacific: 15 storms 
with maximum winds < 65 kt (broken line); 13 typhoons with 
maximum winds < 100 kt (dashed-dotted line); and 14 typhoons 
with maximum winds equal or > 100 kt (solid line). The insert 
contains graphs of the differences in temperatures of the latter 
and the two weaker categories. The rings are illustrated in Fig. 1 . 


6-10, the hurricanes had higher mean T BB . That is, 
the T bb imply that hurricanes have stronger convec- 
tion near the core and stronger subsidence in the 
environment surrounding the storm. Both the con- 
vective and subsidence areas have been observed 
frequently by aircraft reconnaissance and in satel- 
lite imagery (Shenk and Rodgers, 1978). The subsi- 
dence dissipates many of the clouds at distances 
> 650 km from the center and the higher mean 
T bb in these regions are observed because the 
satellite sensor measures the T BB of the sea surface 
in areas of no clouds. 

A similar test was made with the 1970 Pacific 
Ocean tropical cyclones. The composited mean 
T bb for 14 typhoons (current maximum wind V 0 
25 100 kt), 13 typhoons (65 kg «£ V 0 < 100 kt) and 
15 weaker storms (V 0 < 65 kt) are compared in 
Fig. 3. The comparison between the weak storms 
and those with maximum winds > 99 kt is similar 
to that of the 1969 Atlantic storms except for ring 1 
where the T BB for typhoons are higher. This reflects 
the large cloud-free eyes of several of the typhoons 


and the fact that the Nimbus 4 THIR had higher 
resolution than the Nimbus 3 MRIR. If values for 
ring 1 are ignored, intense storms are again colder 
for the inner four rings and warmer at greater radii 
(dashed line. Fig. 3 insert). 

The comparison between the typhoons of moderate 
intensity with the very intense typhoons (solid line, 
Fig. 3 insert), however, gave contradictory informa- 
tion. An examination of the individual cases in- 
volved showed the moderate intensity storms were 
biased toward intensification and the intense storms 
toward weakening or little change. This is especially 
significant because it suggests that the mean T BB are 
also an index of the rate of change of storm intensity . 

The relationship of T BB to rate of change of storm 
intensity was investigated by stratifying the storms 
according to wind intensity changes. Fig. 4 shows 
composited T BB means for cyclones stratified ac- 
cording to intensity changes during the succeeding 
24 hours. Four categories are used: 1) intensifying 
( V Q increase 5= 10 kt), 12 cases; 2) weakening 



Fig. 4. As in Fig. 3 except that the data are stratified accord- 
ing to change of intensity during next 24 h: maximum winds in- 
creasing 10 or more kt (I), maximum winds changing less than 
10 kt (S), maximum winds decreasing 10 or more kt (W), and 
storms which never reached hurricane intensity (T.S.). Only 
storms located south of 30°N were included. The insert shows 
that the intensifying storms have much colder cloud tops within 
8° (888 km) of the storm center than the others. 
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(V 0 decrease s* 10 kt), four cases; 3) little change 
( V 0 change < 10 kt), six cases; and 4) five tropical 
cyclones which did not reach typhoon intensity. The 
cases in categories 1, 2 and 3 all reached typhoon 
intensity. We can note that the storms with the 
greater rate of intensification are associated with the 
lower mean T BB in all rings outward through ring 9. 
The intensifying storms are —18 K colder in rings 
2-4 than the storms that never reached typhoon in- 
tensity, — 10 K colder in rings 1-8 than the weak- 
ening storms and 5- 10 K colder than the storms with 
little change in intensity. 

The time lag between the mean T BB and maximum 
wind of several tropical cyclones examined further 
supports the hypothesis that T BB serves as a predica- 
tor of wind changes. Fig. 5 illustrates a time lag be- 
tween the T bb and the change of the maximum winds 
in 1970 Pacific Typhoons Billie and Hope. The 
lowest mean T BB , which may be said to represent 
maximum convection, of Billie occurred two days 
before the wind maximum was reached and that of 
Hope occurred more than one day prior to the wind 
maximum. Similar data examined in other storms 
suggest that the wind changes lag the T BB changes 
by 24-36 h. 

Fig. 6 is adapted from a simulation experiment 
with a theoretical model by Rosenthal (1978). The 
vertical velocities at 900 mb and the maximum winds 
are plotted against time. Here vertical velocity, 
rather than temperature as in Fig. 5, represents con- 
vection. There is a striking similarity in the time lag 


between maximum convection and the maximum 
winds of the storm in the two illustrations. 

The results of these comparisons of the mean T BB 
and storm intensity strongly suggest that the cloud 
top T bb near the tropical cyclone center may serve 
as an index for both 1) current storm intensity, and 
2) future storm intensity. This dual relationship is 
further explored in the succeeding section where 
equations are developed for predicting storm inten- 
sities 24 h in advance. 

4. Development of predictive equations 

Results of the analyses reported in Section 3 
demonstrate that both current and future storm in- 
tensities are a function of the mean T BB . These re- 
sults also suggest that the relationship between the 
mean T BB and the future intensity varies with at least 
the latitude of the storm, season of the year, sea 
temperatures, lapse rate in the ambient atmosphere 
and past changes in storm intensity. Therefore, it 
appears there is a need to either stratify the data into 
several classes or use several predictors to account 
for the various effects. 

The limited number of cases for which good data 
were available emphasized the need, however, for 
special care in selection of techniques used in devel- 
oping the predictive equations. Instead of using a 
strictly statistical approach in seeking a solution, 
potential predictors were preselected after use of 
a combination of physical reasoning based on re- 



Fig. 5. Temporal changes of mean equivalent blackbody temperatures (rings 
2-5, Fig. 1) and maximum winds of Typhoons Hope (September 1970) and 
Billie (August 1970). The maximum wind changes lag the temperature changes in 
both cases. The temperature scale is inverted to show the lower temperatures at 
the top of the graph. 
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Fig. 6. Comparison of vertical motion and central pressure in a model hurricane showing 
lag with time (Rosenthal, 1978). 


suits of research on tropical cyclones of the last 20 
years and statistical type analyses such as those 
presented in Section 3. Furthermore, because of the 
concern that the use of a large number of predictors 
in the equations would lead to unstable solutions 
which would provide poor forecasts when the pre- 
dictive equations were applied to independent data, 
the arbitrary decision was made that not more than 
four predictors would be used. 4 This required 
special care in selecting the predictors. The mean 
T bb parameter to be used in each case was identified 
by examining data arranged in a systematic fashion 
(e.g., such as in Fig. 4) or by study of previously 
calculated correlation coefficients. 

It was believed that by using only predictors that 
seem to have a close physical relationship with the 
parameter to be predicted that the predictive equa- 
tions developed with even a relatively few cases 
would likely be stable and apply to other storms. To 
further insure that the predictors selected on the 
basis of physical reasoning and statistical deductions 
were truly related with the parameter to be predicted 
and not just correlated by chance, all forecast equa- 
tions developed were tested with independent data 
(i.e., data not used in the preliminary analyses to 
preselect the predictors nor used in the development 
of the regression equations). 


4 In equations tested and used in this report, three predictors 
at most were used. 


Besides the mean T BB parameter, the current max- 
imum wind ( V 0 ) and the change in maximum winds 
during the preceding 24 h ( A V_ 24 ) were used as pre- 
dictors. In addition, parameters putting additional 
emphasis on the degree of organization of the con- 
vection were considered in some of the equations. 
This was usually a parameter measuring the vari- 
ability of the mean T BB in one or more of the rings 
where the subscripts refer to the identifiers of 
rings). Once the predictors had been preselected, 
a modified screening regression approach was used 
in developing the equations in order to determine 
which of the terms contributed most and also to 
determine just how much reduction in variance oc- 
curred as each new predictor was added to the 
equation. 

Stratification of the storms into two groups (V 0 
^ 65 kt and V 0 > 65 kt) produced better results than 
single grouping of all storms. Fifty-eight cases from 
the 1970 western North Pacific tropical cyclones 
were used as dependent data for developing the fol- 
lowing equations. The first three regression equa- 
tions are: 

Weak storms AV+ 24u , = 143.75 - 0.594 

+ 0.389AV-2 4 (2) 

Strong storms V+ 24< = 227.86 - 0.76f lf2 , 3 

+ 0.499AV_ 24 4- 0.398 V 0 (3) 
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All storms V+ Ua = 146.6 - 1.669<r,, 2 

+ 0.855V* - 0.513?,. (4) 

Because of the strength of the relationship of V +24 
with the satellite measured mean TbBj several equa- 
tions were developed for use with satellite data 
alone: 

Strong storms V +24< = 378.51 - 1.225 j li2 , 3 (5) 

Strong storms % V +24 , = 390.72 - 1.246f li2 , 3 

- 0.5060-3 (6) 



y+*« 

AV + „ V 0 

av_ I4 

av_ m 

V « 

-0.547 (-0.472) 
0.504 (0.482) 
0.625 (0.465) 

-0.647 (-0.300) -0.068 (-0.273) 

0.450 (0.213) 

0.177 (0.275) 

-0.291 (0.393) 
0.304 (0.398) 

Reduction in variance of AV +lv dependent data 


Predictors 

Reduction in variance 

7,,; AV_ M 

0.42 0.49 


Table 2. Correlation coefficients for strong storms 
dependent data (independent data). 



v\» 

av +m 

v. 



-0.781 (-0.725) 

-0.589 (-0.642) 

-0.392 (-0.274) 

-0.489 (-0.548) 


-0.056 

-0.286 

-0.197 

-0.407 

AV. t * 

0.662 (0.697) 

0.655 (0.446) 

0.130(0.475) 


v. 

0.551 (0.556) 

-0.191 (-0.144) 


0.130(0.475) 


Reduction in variance of V 

+ 14 , dependent data 


Predictors 


f i^j 

AV_, 4 

tm, AV_ M ; Vo 

Reduction in variance 

0.61 

0.71 

0.77 


5. Discussion of results 

In the weak tropical cyclone set, Table 3 reveals 
that Eqs. (2) and (7) produce better forecast results 
compared to those from persistence techniques not 
only for the dependent data which would be ex- 
pected but also for the independent data sets. The 
difference between regression and persistence (NC) 


Weak storms 

AV +2v = 167.16 

- 0.682 ? 2 ,3 

(7) 

All storms 

V +24a = 200.51 

- 2.213 o- 1>2 

- 0.381 ? 2>3 . 

(8) 


Here the symbols are defined as follows: 

AV +24ip predicted change in maximum wind speed 
(kt) during a 24 h period after satellite 
observation for storm whose current 
maximum winds are ^65 kt 

V + 24 , predicted maximum wind speed (kt) 24 h 

after satellite observation for storm 
whose current maximum winds are 
>65 kt 

V +24a predicted maximum wind speed (kt) 24 h 
after satellite observation for storms of 
any intensity 

T x mean T flfl (K) for area 0- 1 1 1 km about the 
storm center (ring 1) 

t 2f 3 mean T M (K) for area 111-333 km about 

the storm center (rings 2 and 3) 

T i. 2.3 mean T BB ( K) for area 0-333 km about the 

storm center (rings 1, 2 and 3) 

AV _ 24 change in maximum wind speeds (kt) of 

the storm during the 24 h preceding 
satellite observation 

V 0 the current maximum wind speed (kt) in 
the storm at the time of satellite 
observation 

cr 1>2 the standard deviation (K) of the mean 
T bb of all sectors contained in rings 
1 and 2 

cr 3 the standard deviation (K) of the mean 
T bb of all sectors contained in ring 3. 

Table 1 . Correlation coefficients for weak storm cases 
dependent data (independent data). 


Table 1 lists the correlations between the various 
parameters used as predictors with the quantities to 
be predicted, and with each other for the weaker 
storms. Table 2 lists the correlation data for intense 
storms. Eqs. (2) and (7) were developed to predict 
the change in maximum wind speeds during the next 
24 h rather than the maximum wind speed for two 
reasons: 1) the correlation between f 2t3 was higher 
with AV +24 than with V +24 and 2) in the typhoon 
service, operational procedures are such that for 
weak storms the forecaster is believed in many cases 
to have a better estimate of whether a storm is in- 
tensifying than of the absolute value of the maximum 
winds. Of course, by adding V 0 to both sides of the 
equations, one can put them in the same format as 
the other equations. 

The great value of Eqs. (4)-(8) is that no past his- 
tory of the storm is needed. This provides a means 
of making reliable forecasts over isolated ocean 
basins. The current satellite imagery can be used for 
a general classification (over/under 65 kt) and the in- 
frared data may be used to develop the necessary 
T bb parameters. By using the Dvorak technique 
(Dvorak, 1975), one can also obtain V 0 from satellite 
data [Eq. (4)]. 

Errors of forecasts made with the equations were 
compared with those obtained by techniques fre- 
quently used in tropical cyclone intensity forecast- 
ing: persistence (NC) [assumes no change during the 
forecast period], and persistence of change (P) [as- 
sumes the change during forecast period is the same 
rate as the change during the preceding period]. Re- 
sults of this comparison for both the dependent and 
independent data series are summarized in Table 3. 
The independent series were from the 1973 and 1974 
western North Pacific tropical cyclones. 
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Table 3. Mean errors of forecasts. 


Equation 


Dependent data (1970 storms) 


Independent data (1973 and 1974 storms) 

n 

E 

E - NC 

E - P 

n 

E 

E - NC 

E - P 




Weak tropical cyclones (V 0 

^ 65 kt) 




2 

24 

8.6 

-4.0 

-4.2 

33 

12.6 

—4.2* 

-0.6 

4 

24 

10.2 

-2.4 

-2.6 

33 

14.2 

-2.4 

1.0 

7 

24 

9.6 

-3.0 

-3.2 

33 

12.8 

-4.0* 

-0.4 

8 

24 

24.9 

12.3 

12.1 

33 

25.0 

+8.2 

+ 11.8 




Intense tropical cyclones (V fl 

( > 65 kt) 


9 


3 

34 

9.4 

-7.5 

-2.6 

20 

16.1 

-3.1 

-6.9* 

3 (-11 kt) 





20 

11.8 

-7.4* 

-11.2* 

4 

34 

12.1 

-4.8 

0.1 

20 

14.7 

-4.5* 

-8.3** 

5 

34 

11.4 

-5.5 

-0.6 

20 

18.8 

2.0 

5.6 

5 (-11 kt) 





20 

15.6 

-3.6 

-7.4* 

6 

34 

10.8 

-6.1 

-1.2 

20 

18.7 

1.9 

5.5 

6 (-11 kt) 





20 

14.9 

-4.3 

-8.1* 

8 

34 

21.1 

4.2 

9.1 

20 

17.4 

-1.8 

-5.6 




All storms 





2 & 3 

58 

9.1 

-6.0 

-3.2 

53 

13.9 

-3.8* 

-3.0** 

4 

58 

11.4 

-3.7 

-0.9 

53 

14.4 

-3.3* 

-2.5 

8 

58 

22.7 

7.6 

10.4 

53 

22.1 

4.4 

+5.2 


Notes: E is mean error of forecasts made by indicated regression equation. 

E - NC is mean difference in errors of forecasts made by regression equation and by assuming no change in wind speed 
during forecast period. 

E - P is mean difference in errors of forecasts made by regression equation and by assuming the change in wind speed in 
the next 24 h would be the same as the change during the previous 24 h. 

* Differences significant at 1% level. 

** Differences significant at 5% level. 


results in the independent data set is significant at 
the \% level. Eq. (8) yielded only fair results; Eq. 

(4) , however, produced strikingly good results con- 
sidering that satellite data alone were used in fore- 
cast parameters 5 and that this equation applied to 
both weak and strong storms. 

The intense tropical cyclone [Eqs. (3), (5) and (6)] 
provide good forecasts in the dependent data set but 
only Eq. (3) gives results better than persistence with 
the independent data set. However, a bias (see Sec- 
tion 6) is believed to exist in the wind information 
of the independent set and causes the intensity fore- 
casts by the regression equations to be about 1 1 kt 
too high in the mean. When this bias is removed from 
the forecasts for the independent data set, Eqs. (3), 

(5) and (6) are better than persistence. Significance 
at the 1% level is prominent in several compari- 
sons. Eq. (4) shows superiority even without the bias 
removed. 

For all storms combined Eq. (4) gives better re- 
sults than persistence techniques. They are only 
slightly weaker than results from Eqs. (2) and (3) 
combined. This is particularly significant in that fore- 


5 Assumes that V 0 can be obtained by the Dvorak or some 
similar technique. 


casts from Eq. (4) can be made using only satellite 
data. Eq. (8) does not produce results of comparable 
quality. 

The regression equations for the intense and weak 
tropical cyclones have yielded results throughout 
the dependent data set which are -0. 1 to 7.5 kt better 
than the persistence results [Eq. (8) excepted]. The 
independent data set results are -5.6 to 8.3 kt better. 
When the bias error is applied for Eqs. (3), (5) and 

(6) in the independent data set the regression results 
are 3.6-11.2 kt better than the persistence results. 

One of the objectives of this investigation was to 
determine whether the mean T BB were strongly re- 
lated to the current and future wind speeds and 
whether they had predictive value. The fact that the 
regression equations outperform persistence for the 
independent series suggests that the T BB do contain 
predictive information over and above that con- 
tained in the current and previous wind data. Never- 
theless, the fact that equations using AV_ 24 and V 0 
produced superior forecasts to those which used 
only the mean T BB and its derivatives, makes it per- 
tinent to ask if the skill shown in Eqs. (2), (3) and (4) 
is primarily due to the use of wind information. The 
answer to this question is contained in Tables 1 and 2. 

For the weak storms in the dependent data series, 
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the highest correlation with the winds was between 
f 2 , 3 and A V+ 24 . For the independent data the highest 
correlations was with V+ 24 and there were relatively 
little difference in correlations between it and with 
I2.3, AV_ 24 or V 0 . The reduction in variance was 
0.42 when f 2>3 was the sole predictor. Adding A V_ 24 
as a predictor increased the reduction in variance 
to 0.49. 

T 2 ,3 is poorly correlated with AV_ 24 in the weak 
series so they should contribute independently to 
the skill of the regression equation. Note that the 
sign of this correlation changes between dependent 
and independent data. 

For the strong storms dependent data (Table 2) the 
correlation between f, f2(3 with V+ 24 (-0.781) is 
stronger than that between V +24 and either AV_ 24 
(0.662) or V 0 (0.551). While the correlation between 
7 1i2>3 and either V 0 or AV_ 24 is considerably higher 
than for the weak storms, these correlations are 
still lower than the correlations between r lt2t3 and 
either V+ 24 or A V+ 24 (-0.781 and -0.589 vs -0.392 
and -0.489) data. These relationships also hold for 
the independent data. 

In the tests of significance of the differences in 
the errors by the regression equations and the er- 
rors by persistence techniques, consideration was 
given to the size of the samples of the independent 
data. There were 33 cases for the weak category 
selected from 20 different storms from two years. 
There were 20 cases for the strong category selected 
from 9 different storms from two years. Most of the 
cases from the same storm were spaced at intervals 
2* 24 h. There were, however, five cases from one 
storm and three cases from another storm in the 
strong sample which were approximately at 12 h 
intervals. For the weak storms there were five 
storms that contributed two cases and one storm that 
contributed three cases where the interval was ap- 
proximately 12 h. Examination of these cases re- 
vealed that the standard deviation of the errors of 
the regression forecasts were much larger for the 
strong storm cases with the 12 h intervals than for 
the entire sample. For the weaker storms there was 
little difference between the standard deviation of 
the errors for the cases of 12 h intervals and for the 
entire sample. That is, the data from the individual 
cases were sufficiently convincing that there was 
relatively little autocorrelation between the storm 
cases spaced at 12 h intervals and their inclusion did 
not alter the effective sample size. 

6. Bias in forecast results 

The forecast results from Eq. (3) are better for 
the independent data than those obtained from per- 
sistence. However, examination of the errors re- 
veals that the regression equation forecasts maxi- 
mum winds 11.3 kt too high in the mean for the 


independent cases of 1973 and 1974. This bias was 
calculated by taking the algebraic mean 6 of the fore- 
cast errors. The equations and data were tested to 
determine the cause of the bias. The following were 
checked for possible biasing effects: 1) satellite sen- 
sor calibration variation; 2) time of year of storm 
occurrence; 3) location of storm; and 4) operation 
and best track estimates of maximum wind speeds. 

To check the sensitivity of results from the equa- 
tions to possible bias corrections, the constant term 
of Eq. (3) was reduced by various amounts ranging 
from 8 to 17 kt. This variation caused only minor 
variations in the mean of the absolute values of the 
errors for the independent series (12.7- 1 1 .6 kt which 
compares with 16.1 kt when no bias was applied). 
This suggests that a bias correction in the range 8- 17 
kt that was well substantiated might be used to adjust 
the constant in Eq. (3) when it was used with data 
from the 1973 and 1974 seasons. Other equations 
were tested with similar results. 

The biasing effect caused by using T BB observa- 
tions from different satellites was minimum. The 
upper and lower limits of the mean T BB values and 
their range in the tropical cyclones as measured in the 
storms by Nimbus 3, 4 and 5 agreed within <^2°C. 

Seasonal, latitudinal and longitudinal variation be- 
tween the dependent and independent data sets re- 
sulted in differences which would account for only a 
few knots (-2 to 4 kt) difference in forecast for the 
dependent and independent data and even then only 
for a few storms. The effect of these factors is not 
believed, therefore, to be of sufficient magnitude to 
account for the large bias found to exist. 

The operational estimates of maximum wind, 
however, could account for the bias. Atkinson and 
Holliday (1977) discuss the “considerable uncer- 
tainty involved in existing equations” for estimating 
maximum winds and report changes in procedures 
at the Joint Typhoon Warning Center at Guam 
(JTWC) that may affect results of our experiment. 
Table 4 compares the maximum wind calculated for 
a range of central pressures ( P c ) from the following 
three equations (Atkinson and Holliday, 1977): 

Equation by Atkinson and Holliday 

V mi = 6.7(1010 - P c )°- 644 (9) 

Equation by Takahashi 

V*. = 13.4(1010 -P f ) 05 (10) 

Equation by Takahashi (high lat.) 

Vm, = 11.5(1010 ~P c y \ (11) 

Here P c is the central pressure (mb) of the storm 
and V m the maximum wind speed (kt). When these 

* The errors listed in Table 3 are all means of absolute values of 
individual errors. 
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Table 4. Relations between minimum central pressure and 
maximum winds in tropical cyclones. 


Pc 

v mt 

V- 

V„ 

1000 

30 

42 

36 

990 

46 

60 

51 

980 

60 

73 

63 

970 

72 

85 

73 

960 

83 

95 

81 

950 

94 

104 

89 

940 

103 

112 

96 

930 

113 

120 

103 

920 

122 

127 

109 

910 

130 

134 

115 

900 

138 

141 

121 

890 

146 

147 



equations are applied to the independent data set 
cases the mean difference between the results from 
Takahashi’s and Atkinson and Holliday’s equations 
is near 10 kt. Holliday 7 reports that Eq. (9) was used 
for developing most of the “best track” informa- 
tion for the 1973 and 1974 years while Eqs. (10) and 
(11) were used in 1970. These results suggest that 
the changes in procedures at JTWC for estimating 
wind maxima may account for most of the bias found. 

To further test the hypothesis that the bias is due 
to change in procedures for obtaining maximum 
winds, the 1973, 1974 data were used dependently 
to develop equations similar to Eqs. (3) and (6). 
These new equations 

V +24f = 243.91 - 0.8857^2,3 

+ 0.385 A V_24 + 0.424 V 0 > (12) 

V +24 = 395.3 - 1.333 - 0.404,3, (13) 

were substracted from Eq. (3) and (6), respectively, 
and the results simplified to provide 

V +24 , - V +24j 

= 8.5 + 0.125(7 1 , 2j3 - 210) 

+ 0.1 14A V_ 24 - 0.026( V 0 - 65), (14) 

V+ 24 , - V +24 , = 13.078 + 0.087(7^2,3 - 210) 

- 102(o-3 - 6). (15) 

When Eq. (14) is solved using representative values 
(7 1 ,2,3 = 230 K, V 0 = 90 kt) the result is 10.9 kt. 
When Eq. (15) is solved using representative values 
(7 1>2 ,3 = 230 K, cr 3 = 10 K) the result is 13.7 kt. 
Forecasts made by Eq. (3) with the constants re- 
duced by 11.3 and those made with Eq. (12) have a 
correlation coefficient of 0.99. These various num- 
bers compare favorably with the 11.3 kt bias found 
and are all within the range 8-17 kt. Thus, the re- 


sults lend strong support to the contention that much 
of the bias in Eq. (3) exist due to the change in pro- 
cedures for estimating wind maxima when preparing 
the “best track” information at JTWC. 

7. Simulated operational test results 

Simulated operational tests were made using the 
regression equations (2), (3) and (4). Wind data from 
typhoon advisories for the 1973, 1974 storms were 
incorporated as input parameters in place of best 
track winds. The results produced even greater in- 
creases in errors in persistence forecasts than in the 
regression forecasts. The regression results could 
not be readily compared with the JTWC forecasts 
for these test cases. JTWC advisory (forecast) times 
did not coincide exactly with satellite observation 
times and not all JTWC forecasts were available. 
To the extent information was available, it was pos- 
sible to determine that the JTWC forecast errors of 
the test cases greatly exceeded their mean absolute 
error of 13.5 kt for the entire 1973 and 1974 seasons. 
The JTWC forecast errors were also larger than the 
mean errors of the regression equations for the cases 
where adequate comparisons could be made. It ap- 
pears from these results as well as those discussed 
in Section 5 that the regression equations will yield 
better forecasts than other techniques being used. 

8. Conclusions 

This investigation demonstrated that there is a 
strong relationship between satellite measured T BB 
of cloud tops near the core of the tropical cyclones 
and both the current and future storm intensity. This 
was demonstrated for four different seasons of data 
(one in the Atlantic and three in the Pacific area). 
The statistical correlation, moreover, reveals that 
the future intensity (V +24 ) is more closely related 
than the current intensity (V 0 ) to the mean T BB 
(correlation -0.781 vs 0.392, Table 2). 

The regression equations provided forecasts 
superior to those of persistence techniques when 
applied to independent data. The equations pro- 
duced better results than other techniques in simu- 
lated operational tests. The results are sufficiently 
convincing to warrant use of the equations as ex- 
perimental forecast tools. It is particularly important 
that such good results are obtained from quantified 
T bb data, a parameter readily obtainable from a 
single satellite observation. 

Further testing is recommended before imple- 
menting the technique in operational forecasting 
routines because of the limited sample size which 
was available in this investigation. It would be 
desirable to expand the investigation to include more 
cases as well as other tropical cyclone basins. The 
data contained in this report appear convincing, 


7 Personal communication. 
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however, that there is a relationship between T BB 
and future maximum winds in tropical cyclones that 
is strong enough to have predictive value. 
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